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Abstract. We develop a Hamiltonian theory of a time dispersive and dis- 
sipative inhomogeneous medium, as described by a linear response equation 
respecting causality and power dissipation. The proposed Hamiltonian couples 
the given system to auxiliary fields, in the universal form of a so-called canon- 
ical heat bath. After integrating out the heat bath the original dissipative 
evolution is exactly reproduced. Furthermore, we show that the dynamics as- 
sociated to a minimal Hamiltonian are essentially unique, up to a natural class 
of isomorphisms. Using this formalism, we obtain closed form expressions for 
the energy density, energy flux, momentum density, and stress tensor involving 
the auxiliary fields, from which we derive an approximate, "Brillouin-type," 
formula for the time averaged energy density and stress tensor associated to 
an almost mono-chromatic wave. 



1. Introduction 

The need for a Hamiltonian description of dissipative systems has long been 
known. Forty years ago Morse and Feshbach gave an example of an artificial Hamil- 
tonian for a damped oscillator based on a "mirror-image" trick, incorporating a 
second oscillator with negative friction [30]. Ch 3.2]. The resulting Hamiltonian is 
tm-physical: it is unbounded from below and under time reversal the oscillator is 
transformed into its "mirror-image." The artificial nature of this construction was 
described in [301 Ch. 3.2]: "By this arbitrary trick we are able to handle dissipa- 
tive systems as though they were conservative. This is not very satisfactory if an 
alternate method of solution is known..." 

We propose here a quite general "satisfactory solution" to the general problem 
posed in [30] by constructing a Hamiltonian for a time dispersive and dissipative 
(TDD) dynamical system without introducing negative friction and, in particular, 
without "mirror-images." Developing a Hamiltonian structure for a TDD system 
might seem a paradoxical goal — after all, neither dissipation nor time dispersion 
occur in Hamiltonian systems. However, we will see that if dissipation is introduced 
via a friction function, or susceptibility, obeying a power dissipation condition — 
as it is for a linear dielectric medium described by the Maxwell equations with fre- 
quency dependent material relations — then the dynamics are exactly reproduced 
by a particular coupling of the TDD system to an effective model for the normal 
modes of the underlying medium as independent oscillating strings. For the com- 
bined system we give a non-negative Hamiltonian with a transparent interpretation 
as the system energy. 
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An important motivation behind this effort is the clarification of the definition 
of the radiation energy density and stress tensor for a dissipative medium in the 
linear response theory, e.g., a dielectric medium with complex valued frequency 
dependent material relations. An intrinsic ambiguity in this definition has led to 
problems interpreting the energy balance equation [22], Sect. 77], [6] Sect. 1.5a], 
[Ell Sect. 6.8], [26]. These difficulties do not persist if a fundamental microscopic 
theory is considered. Consequently a number of efforts [26 ] I23 | 132] have been made 
to construct a macroscopic theory of dielectric media, accounting for dispersion and 
dissipation, based on a more fundamental microscopic theory. It might seem that 
the introduction of an explicit realistic material medium is the only way to model 
a TDD system. However, the construction of this paper shows this is not so and 
provides a consistent macroscopic approach within linear response theory. 

As an example of our general construction, we analyze here TDD dielectric me- 
dia, including a detailed analysis of the electromagnetic energy and momentum 
densities. Part of that analysis is the derivation of an approximate formula for the 
time averaged Maxwell stress tensor similar to the Brillouin formula for the time 
averaged energy density [22l Section 80] . 

Another important benefit of the approach developed here — and in our previous 
work [7] — is that the present formulation allows to apply the well developed 
scattering theory for conservative systems |33j to the long long standing problem 
of scattering from a lossy non spherical scatter — analyzed by other methods with 
limited success [29] . This will be discussed in detail in forthcoming work [9] . 

1.1. Dissipative systems. We consider a system to be dissipative if energy tends 
to decrease under its evolution. It is common, taking energy conservation as a 
fundamental principle, to view a dissipative system as coupled to a heat bath so 
that energy lost to dissipation is viewed as having been absorbed by the heat bath. 

We have shown in [7] that, indeed, a general linear causal TDD system can 
be represented as a subsystem of a conservative system, with the minimal such 
extension unique up to isomorphism. Let us summarize the main ideas of that 
work here. One begins with an evolution equation accounting for dispersion and 
dissipation, of the form 

(1.1) dtf (t) = Lf (t) - ^ a (r) / (t - t) dr + r(t), 

Jo 

where / (t) describes the state of the system at time t, specified by a point in a 
complex Hilbert space Hq, and 

(1) L = — L4 with A a self-adjoint operator on Hq, 

(2) a (t) is an operator valued function, called the friction function and as- 
sumed to be of the form 

(1.2) a (t) = aoo 6 (r) + a (r) , 

wtih a (t) strongly continuous for r > and self-adjoint. 

(3) r(t) is an external driving force. 

For a(t) satisfying a power dissipation condition [ ()1.36j) . below] one then constructs 
a complex Hilbert space Ti., an isometric injection / : Hq — > TLq and a self adjoint 
operator A on Tt such that the solution fit) to (|1 . 1[) equals the projection V F{t) 
onto Hq of the solution F(t) to 

(1.3) d t F(t)=£F{t) + 7r(t), 
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with C = — L4. That is 

(1.4) /(<) = PF(t) = f J t e (t - t ' )£ /r(t')di'. 

J — oo 

The main object of this work is to extend [7] by considering dissipative systems 
with a given Hamiltonian structure. Here and below we call a system Hamiltonian if 
its phase space is endowed with a symplectic structure and a Hamiltonian function 
such that: i.) the system evolves by Hamilton's equations, and ii.) the physical 
energy of the system in a configuration associated to a phase space point u is equal 
to the value of the Hamiltonian function at u. Accordingly, a dissipative system 
is by definition not Hamiltonian. Nonetheless, almost every dissipative system of 
interest to physics is a perturbation of a Hamiltonian system, with the perturbation 
accounting for dispersion and dissipation. 

The property of being Hamiltonian, as defined above, is more than a formal 
property of the evolution equations, as it also involves a physical restriction equating 
the Hamilton function and the system energy. For a linear system such as (|1.3[) 
there are many ways to represent the evolution equations as Hamilton's equations. 
We circumvent this ambiguity by supposing given a Hamiltonian structure on the 
given TDD system, whose evolution is a suitable perturbation of the Hamilton 
equations. We then ask, and answer affirmatively, the question, "Is there a natural 
way to extend the given Hamiltonian structure to the unique minimal conservative 
extension of [7] so that the extended system is Hamiltonian?" This will be achieved 
in a self contained way below, with reference to the extension of [7] , by constructing 
a Hamiltonian extension with additional degrees of freedom in the universal form 
of a canonical heat bath as defined in [14j Section 2], [37, Section 2]. 



1.2. Hamiltonian systems. We suppose given a dynamical system described by 
a coordinate u taking values in phase space, a real Hilbert space V. On V there is 
defined a symplectic form 3(u, v) — (u, Jv), with J : V — > V a linear map such that 

(1.5) J T J = 1, J 2 = -l. 

We call a map J satisfying (|1.5p a symplectic operator. Throughout we work with 
real Hilbert spaces and use M T to denote the transpose of an operator M, i.e., the 
adjoint with respect to a real inner product. Additional notation, used without 
comment below, is summarized in Appendix along with the spectral theory for 
operators in real Hilbert spaces. 

The evolution equation (in the limit of zero dissipation) is to be Hamiltonian 
with respect to 3- Thus, we suppose given a Hamiltonian function h(it) such that 
when dissipation is negligible u evolves according to the symplectic gradient of h 

(1.6) d t u = J— — . 

ou 

For most applications the Hamiltonian h (u) is the system energy and is nonnegative 
(or at least bounded from below). Mostly, we consider a quadratic, non- negative 
Hamiltonian 



(1.7) 



h (u) = i (Ku , Ku) 
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leading to a linear evolution equation. However, there is a natural extension of the 
results presented here to a nonlinear system with 

(1.8) h(u) = -(Ku,Ku)+ti(u), 

where h' (u) is an arbitrary function of u. The construction carries over, provided 
the dissipation enters linearly as in (|1.10b'[) below. To keep the discussion as simple 
as possible — and to avoid the difficult questions of existence and uniqueness for 
non-linear systems — we consider the Hamiltonian (|1.7p throughout the main text. 
A few examples illustrating this point are discussed in Appendix [Cl 

We call the operator K the internal impedance operator (see (j!.10b[) below), and 
suppose it to be a closed, densely defined map 

(1.9) K : V{K) -> H, V(K) C V 

with H the stress space. The (real) Hilbert spaces V and H are respectively the 
system phase space and the state-space of internal "stresses. "Q The space of fi- 
nite energy states is the operator domain T>(K). Physical examples and further 
discussion of the operator K are given in Section [3] Technical assumptions and a 
discussion of the dynamics on T>(K) are given in Section [5] 

The equation of motion, in the absence of dissipation, is obtained from (|1.61 II. 7p 
by formal differentiation. It is convenient split the equation in two: 

(1.10a) dtu{t) — JK T f{t) (evolution equation), 

with 

(1.10b) f(t) = Ku{t) (material relation without dispersion or dissipation). 

When dissipation is included, we replace (jl.lObj) with a generalized material rela- 
tion, 

poo 

(1.10b') /(*)+/ dT X (r)f(t-r) = Ku(t), 

Jo 

where \ is the operator valued generalized susceptibility, a function of r > with 
values in the bounded operators on H. Note that the integral in (|1.10b'[) explicitly 
satisfies causality: the left hand side depends only on times t — r < t. 

The structure of (jl.lOal Il.l0b'[) mirrors the Maxwell equations for the electro- 
magnetic (EM) field in a TDD medium. For a static non-dispersive medium — 
see Section 2] — eq. (|1.10al) and (|1.10b[) correspond respectively to the dynami- 
cal Maxwell equations and the material relations. (The static Maxwell equations 
amount to a choice of coordinates.) Dispersion and dissipation are incorporated in 
(|1.10a[ ll.lOb'j) by modifying the material relations in the same fashion as in the 
phenomenological theory of the EM field in a TDD medium. 



Abstractly, it is not strictly necessary to distinguish V and H. We could replace H by V 
and K by \K\ (see Appendix |A]| . However, that might be physically unnatural, and we find that 
the distinction clarifies the role of dissipation and dispersion in applications. In particular, the 
impedance operator is dimensionful (making it necessary to distinguish domain and range) unless 
we parametrize phase space by quantities with units ^/energy. 

From a mathematical standpoint, using \K\ may introduce complications. For instance, 
with V = L 2 (R 3 ;C), H = L 2 (R 3 ;C 3 ) and Ku(f) = Vu(f), the associated Hamiltonian, 
h(u) = J K 3 d 3 r |Vn(r)| 2 , produces the evolution dtut(r) = — iAut(r), taking J = multiplica- 
tion by i. Of course we might take H = V and K = tJ—A = |V|, instead. But it is more elegant 
(and more natural) to work with the differential operator V. 
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The vectors u and / of the TDD system ([1.10a[ ll.lOb')) may interpreted phys- 
ically as follows: u specifies the state of the system and / specifies internal forces 
driving the dynamics. Thus, we refer to / as the kinematical stress. Similarly, we 
refer to Ku as the mechanical or internal stress, as the square of its magnitude is 
the energy of the system. In a dispersionless system, these quantities are equal, 
but in a TDD system they are not, and are related by (jl.lOb'j) , incorporating time 
dispersion 

Associated to the dispersionless system (jl.lOal II . 10b|) is the initial value problem 
(IVP), which asks for u(t), t > to given the initial condition u(to) = uq. Under 
suitable hypotheses on K and J this problem is well-posed for uo £ T>(K), with 
existence and uniqueness of solutions provable by standard spectral theory (see 
£15. 1|) . However, for the TDD system (|1.10al[T.10b'|) . the initial value problem is not 
well defined, because the integral on the l.h.s. of (|1.10b'|) involves f(t) for t — ► — oo. 
This dependence on history forces us to ask, "how were the initial conditions /o and 
wo produced?" Thus a more physically sound approach is to suppose the system is 
driven by a time dependent external force p(t) (which we controll). This leads us 
to the driven system: 

(1.11a) d t u{t) = JK T f(t) + p(t) 

(1.11b) Ku(t) = f(t)+ / dT X (r)f(t-T), 

Jo 

with initial conditions 

(1.11c) lim u(t)ekeiK, lim f(t) = 0, 

t — > — oo t — •> — oo 

so at t = — oo the system was at rest with zero energy. In the absence of dispersion, 
when x = 0, eqs. (jl.llal ll.llbj) reduce to 

(1.12) d t u{t) = JK T Ku{t) + p{t). 

It is useful to note that (|1.12|) is Hamilton's equation for the time dependent Hamil- 
tonian h t (u) = h(it) — (Jp(t),u). 

We shall generally take the external force to be a bounded compactly supported 
function p : K — > V. More generally we might ask only that p £ L 1 (M, V) or 
even allow p to be a measure. The initial value problem for (|1.12p amounts to the 
idealization pit) — uo5{t — to). 

1.3. Hamiltonian extensions. The main question addressed here is: when does 
the system described by ([1.11)1 admit a Hamiltonian extension? We restrict our- 
selves to looking for a quadratic Hamiltonian extension (QHE), defined below. Our 
main result is the existence of a QHE under physically natural conditions on the 
susceptibility: 

Theorem 1.1. Under mild regularity assumptions for the system operators K and 
X (spelled out in Section [5]), if x is symmetric, 

(1.13) X (t) T = X(t) , 

We could consider a relation inverse to expressing the kinematical stress as a function 

of the mechanical stress, f(t) = Ku(t ) + J n °° dr x(r) Ku(t — r). Under the power dissipation 
condition, Q1.140 below, we may invert l|l.l0b'[l to obtain this equation and vice versa. However 
appears in the standard form of Maxwell's equations and is most convenient for our 

analysis. 
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then there exists a quadratic Hamiltonian extension of the system if and only 

if x satisfies the power dissipation condition (PDC) 

1 



Im {C£(C)} = 57 {Cx(0 ~ CxiCY }>OforallC = u+ir l> r ] > 0, 



(1.14) 

with x the Fourier- Laplace transform of x, 



(1.15) 



dte^ X (t). 



Remarks: i.) The operator x(£) is complex linear, defined on the complcxification 
CH of the real Hilbert space H (see Appendix [A"|) . As indicated, the imaginary 
part in (| 1 . 14[) refers to the imaginary part with respect to the Hermitian structure 
on CH. Due to the symmetry condition (|1.13p . this is also the imaginary part with 
respect to the complex structure, i.e., ImcCx(C) = h {CxiC) ~ C*x(C)*}i where 
* denotes complex conjugation, A*v — (Av*)* . ii.) As mentioned above, the 
result extends with no extra effort to a non-linear system, with a non-quadratic 
Hamiltonian h(u), provided the dissipation is introduced linearly as in (|1.10b f |) . In 
that case, the extended Hamiltonian is, of course, not quadratic as it maintains the 
non-quadratic part of the initial Hamiltonian h(u). 

We verify the theorem by constructing an explicit extension based on the follow- 
ing operator valued coupling function 



duj e 



(1.16) 

and the associated map 

(1.17) Tip = 

The extended Hamiltonian is 
1 



Z7T 



dm cos(ws) y 2wlmx(uj) , 



dos(o-)ip(a) , T : i 2 (R, H) -> H. 



(1.18) 

with 
(1.19) 



E(U) = 



\Ku-T(p\\ H 



\ds<fi(s)\\ H 



~{KU,KU), 



KU = 




The extended impedance K, is a densely defined closed map from extended phase 
space 

(1.20) V = V (B L 2 (R,H)(B L 2 {R,H) 
into extended stress space 

(1.21) W := H(BL 2 (R,H)®L 2 (R,H). 

The symplectic structure on V is given by the following extension of J: 

/ J 

(1.22) J = lo 1 | V. 

loi o 
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We denote by Iy and Ih the isometric injections V V and H <—> H respectively: 

(1.23) I v u = I Oj and I H f = 

The driven Hamilton equations for the extended system are 

(1.24) d t u{t) = JK T f(t)+p{t) 

(1.25) d t 6(s,t) = d 2 Ms,t)+i(s)f(t), 

(1.26) dMs,t) = 9{s,t) 

where we have set the driving force R(t) = Iyp(t) and introduced the kinematical 
stress / in terms terms of the extended system: 

/oo 
da<;(<j)<p{<j,t) . 
-oo 

We think of ip as the displacement of an infinite "hidden string" in K x H . The equi- 
librium configuration of this string is 1 x {0}, and displacements in the directions 
described by H move harmonically driven by the time dependent force ?(s)/(t). 

This explicit extension is an example of what we call a Quadratic Hamiltonian 
extension of . Namely, it is a dynamical system described by a vector coordi- 

nate U, taking values in an extended phase space V, with the following properties: 

(1) The system is a quadratic Hamiltonian system. That is, there are an ex- 
tended symplectic operator J : V — > V and an extended impedance op- 
erator JC : V — > H, taking values in extended stress space, such that the 
evolution of U £ V is governed by 

(1.28a) d t U{t) = JIC T F(t) + R(i), 

(1.28b) F{t) = fCU(t) 

with R(t) the external force. In other words, the dynamics are Hamiltonian 
with symplectic form 2{U, U') = (U^U 1 ) and Hamiltonian 

(1.29) H(t/) = ^(KU,KU). 

(2) The system extends (|1.11[) , as follows. There are isometric injections 

(1.30) Iy-.V^V and I H :H^H 
such that 

(1.31) I H K = KI V , JI V = IvJ, 

and the solution u(t) to (jl.lip . with given initial condition w_oo £ ker K 
and driving force p(t), is u(t) = IyU(t), where U(t) is solves (|1.28p with 

(1.32) lim U(t) = Iyu-co, and R(t) = I v p(t). 

t — > — oo 

Thus the TDD dynamics of p. lip may be modeled by describing u(t) as one com- 
ponent of an extended vector. The motion of the extended system is reversible, but 
an irreversible motion of the underlying TDD system results. This is demonstrated 
in its simplest form by the Lamb model [19] — see Fig. Q] — in which the energy 
of an oscillator escapes to infinity along an attached flexible string. For a simple 
damped harmonic oscillator, the Hamiltonian theory proposed here is precisely the 
Lamb model, and is otherwise a generalization of the Lamb model, obtained by 
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Figure 1. The Lamb model introduced in [TH] to describe radia- 
tion damping, is a point mass attached to an infinite elastic string 
and a Hook's law spring. The point mass evolves as a classical 
linearly damped oscillator. 

coupling an infinite classical elastic string to every degree of freedom of the initial 
Hamiltonian system, illustrating that, from the standpoint of thermodynamics, dis- 
sipation in classical linear response is an idealization which assumes infinite heat 
capacity of (hidden) degrees of freedom. 

1.4. Evolution in stress space and a minimal extension. The extension just 
described is closely related to the extension theory of 7 , summarized in ill. II To 
understand the relation between the present work and [7], it is useful to recast the 
evolution (jl.lip in stress space. If x(r) is, say, continuous on [0, oo) and differen- 
tiable for r > then, by [jl.llbj) . 

(1.33) d t Ku(t) = d t f(t)+ / dr X (r)d t f(t-T)dT 

Jo 

= dt f(t) + X (0)f(t) + / dr [d T x} (r)f(t - r)dr . 
Jo 

Combining this with (|l.lla|) . we obtain: 

(1.34) dtf{t) = KJK T f(t) — f dTa{ T )f(t-T)+Kp(t) , 

J[0,oa) 

where a is the operator valued distribution 

(1.35) o(t) = X (0)S(t) + d tX (t) . 

The evolution (|1.34p is essentially of the form with the minor difference 

that it is defined on a real Hilbert space with a skew-symmetric generator. This is 
of no consequence, as the main result of [7] holds in this context: 
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Theorem 1.2. Suppose given a linear dynamical system described by a point f 
taking values in a real Hilbert space H which evolves according to with a 

skew- symmetric generator L — —L T . If the friction function a(-) — ctoo5(-) + a(-) 
satisfies the power dissipation condition 

/OO POO 
/ (f(t),a e (t-T)f(T))dtdT>0, 
-oo J— oo 

for compactly supported continuous functions f : R — > Hq } where 

(1.37) a e (t)=2a 00 5(t) + { , - oo < i < oo, 

then there exist a real Hilbert space extension H «— 5 and a skew- symmetric oper- 
ator C defined on TL such that (|1.4p holds. 

If, furthermore, the pair (7i, C) is minimal, in the sense that 7i is the smallest 
invariant subspace for C containing the range of I , then the pair (Ti., C) is unique 
up to transformation by an orthogonal isometry. 

Remark: The existence of an extension follows from the results of [7] applied to 
the complexification of (|l.ip fwith / a point in <CH). This extension will not be 
minimal in general, but we can restrict the generator to a suitable real subspace 
to get the minimal extension. Uniqueness may be verified by the arguments of [7] . 
For completeness we give a more detailed sketch of the proof in Appendix [Bl 

The power dissipation condition (|1.14p of the present work implies the PDC 
(|1.36[) of [7] for the friction function a defined in Q1.35p . since 

(1.38) a e (t) = d tX °(t), 
with x° the odd extension of the susceptibility %, 

(1.39) X °(t) := 

In the present work, as in [7], the energy of the dissipative system at time t is 
£{t) :— \\\f{t)\\ 2 - For a trajectory /(•) which evolves according to (|1.34|) this gives 
a total change in energy from t = — oo to t = +oo of 

/OO P oo 

d t £(t)dt = Wfr[/]+ / Rc(f(t),r(t))dt, 

where 

i rOG pOO 

(1.41) W fr [/]:=-- / / (f(t),a e (t-T)f(T))dtdT, 

^ J -oo J —oo 

It is natural to interpret the two terms on the r.h.s. of (|1.40[) as the total work done 
by the frictional and external forces, respectively. Thus, the PDC (|1.36p essentially 
requires that the total work done by frictional forces is always non-positive. This 
physically natural property also provides a necessary and sufficient condition for an 
extension of the form (|1.3|) to exist. 

The theorem guarantees the existence of a unique minimal extension of the form 
(|1.3[) to the evolution in stress space (|1.34[) . However, for this to be a Hamiltonian 
extension we must impose a Hamiltonian structure on the dynamical system (| 1 . 3[) . 
In particular we must express the generator £ as a product 

(1.42) C = JCJJC T , 
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with J a symplectic operator. But, given a skew-adjoint operator, there are in 
general many ways to decompose it in this fashion and thus many ways to impose 
a Hamiltonian structure on the evolution (jl.3p . For the resulting structure to be 
naturally related to the Hamiltonian structure of the original dynamical system 
(|1.10all.l0b[) it is necessary that K and J extend the original operators K and J 
respectively. The main point of this work is to exhibit an explicit Hamiltonian ex- 
tension with these properties, that may then be used in the analysis of conservation 
laws for the dissipative system 

So, by following the motion of the extended stress vector F(t) = JCU(t), we find 
one extension of the type guaranteed by Theorem 11.21 Indeed, f| 1 . 28[) implies 

(1.43) d t F(t) = ICJIC T F(t) + I H K P (t), 

where the generator ¥LJK? is skew-symmetric and has the formal expression 

(KJK T 

(1.44) K.JK? = T T 

V o 

The solution to (] 1 .43[) is easily expressed 

r-t 

(1.45) F(t) = / e^- 1 '^^ I H Kp{t')dt', 




in terms of the one parameter group e**" 1 -^ of orthogonal transformations. By the 
properties of the QHE, the solution to (| 1 .34|) is therefore expressed as 

(1.46) f(t) = I T H F(t) = f ll^-^^I H K P {t')At'. 



It is natural to ask is the extension (|1.44|) is the unique minimal extension of 
Theorem 1 1.2 1 In fact, it is not minimal. Indeed, easily shows that any configuration 
of the hidden string resulting from a physical driving force Iyp(t) is symmetric 
under s <-> — s. That is, we still have a QHE if we replace V and 7i respectively by 

(1.47) V s = V © <S(R, H) © <S(M, H) and H s = H © <S(R, H) © A(R, H), 

with 

S(R,H) = {0eL 2 (R,H) : 

A(R,H) = {</>€ L 2 {R,H) : cj)(s) =-</)(- s)}. 



Note that K : V s — > Tt s and that J : V s — > V s ■ If the kernel of the susceptibility 
kerx(w) is non-trivial on a set of positive measure we will see in £12.11 that further 
reductions are possible. 

There is however no harm in working with an extension which is non-minimal, 
which we do for convenience of notation. Indeed, by (| 1 .45[) . the solution F(t) 
remains in the subspace TCq that is the smallest invariant subspace for /Cj7X T con- 
taining the range of Ih- The restriction of KJK7 to this subspace is the unique 
minimal extension of Theorem 11.21 Thus even if we employ a non-minimal exten- 
sion, we effectively work with the unique minimal extension anyway. In ^2.11 we 
give an explicit description of 7io as well as the minimal subspace Vo C V such that 
K : V -» Ha and J : V -» V . 

Finally, we note that even in the minimal extension there is a great deal of free- 
dom to change variables and thus alter the explicit expressions for the extended 
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impedance fC. Indeed, given a symplectic operator J there is a natural symplec- 
tic group of symmetries of phase space V consisting of linear maps M. such that 
Ai T J^Ai = J . Likewise the Hamiltonian i||/C[/|| 2 does not change if we replace 
the impedance JC by OK. with O any orthogonal map of stress space, O t O = 1. 
Thus the impedance K is essentially defined only up to rc-paramctcrizations of the 
form 

(1.49) K ' ^ O/CA/T 1 , O t O = 1 and M T JM = J. 

We refer to a combined mapping (|1.49p of phase and stress space as a symplec- 
tic/ 'orthogonal isomorphism. (A symplectic map Ai T J^Ai = J need not be bounded 
in infinite dimensions, making it somewhat difficult to formulate the change of vari- 
ables (|1.49[) in complete generality.) 

1.5. Relation with the previous literature. Analysis of a dispersive and dis- 
sipative medium based on the construction of its Lagrangian or Hamiltonian is a 
well established area, see [33J 123 1221 EZ1 HE] arid references therein. However, all of 
those works have relied on specifying an underlying micro-structure for the material 
medium, such as an infinite lattice of dipoles as in [23j . In contrast, our approach 
is phenomenological. Our hidden variables are not "real" microscopic variables as 
in [53J , but describe effective modes which exactly produce a prescribed causal fre- 
quency dependent susceptibility. As regards the underlying microscopic theory, our 
construction can be seen as giving an effective Hamiltonian for those modes well 
approximated by linear response. 

In this section, we compare the approach developed in this paper and our previ- 
ous work [7] with a number of other efforts to describe dissipative and or dispersive 
media via extensions (instead of microscopic variables). 

1.5.1. Dilation theory. The dilation theory — beginning with the Sz.-Nagy-Foias 
theory of contractions j43j [44] and Naimark's theory of positive operator valued 
measures |31j and subsequently extended by a number of other authors — was the 
first general method for constructing a spectral theory of dissipative operators and 
has ultimately provided a complete treatment of dissipative linear systems without 
dispersion. A key observation of our previous work [7j is that many of the classical 
tools of dilation theory, in particular Naimark's theorem, are useful for describing 
the generic case of dissipative and dispersive systems. 

Let us recall the basics of the dilation theory as presented by Pavlov in his 
extensive review |34| and his more recent work [35] . Although there are a number 
of approaches to the subject, Pavlov uses Lax-Phillips scattering theory, [25], which 
provides a conceptually useful picture of the extended operators. That theory 
assumes the existence of: (i) a dynamical unitary evolution group Ut = e lf2t in a 
Hilbert space H where f2 is a self-adjoint operator in H; (ii) an incoming subspace 
D- C H invariant under the semi-group Ut, t < 0, and an outgoing subspace 
D + C H invariant under to the semi-group Ut, t > 0. The invariant subspaces 
D± (called scattering channels) are assumed to be orthogonal. Then one introduces 
the observation subspace Q — H Q (D_ ©£)+), assumed to be co-invariant with 
respect to the unitary group, in the sense that the restriction of Ut, t > 0, to Q is a 
semigroup, i.e., 

(1.50) Z t = P Q U t \ Q = e iBt , 
where Pq is the orthogonal projection onto Q. 
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In many interesting cases the generator B of the semigroup Zt is dissipative, 
i.e. ImB > 0, and the relation (|1.50j) provides a natural setting for dissipative 
operators within the Lax-Phillips scattering theory. The dilation theory reverses 
the Lax-Phillips construction by constructing the Lax-Phillips spaces given B and 
Q. The generator fl of the constructed unitary group is called the dilation of B 
and has the property 

(1.51) f(B) = p Q m)\ Q 

for suitable analytic functions /. Thus the self-adjoint operator f2 provides an 
effective spectral theory for the non-self adjoint B. 

Unfortunately, the dilation theory fails to describe many important physical sit- 
uations simply because the assumption that dissipation occurs without dispersion, 
i.e. that Zt is a semi-group, is too restrictive. In systems such as (jl.lj) dissipation 
comes with dispersion, and dilation theory applies only in the very special case of 
instantaneous (Markovian) friction a(t) — a^d (t). Many phenomenological mod- 
els, such as Lorentz or Debeye dielectric media, employ friction functions which 
are not instantaneous. For such systems one must use a more general approach as 
developed in [7] and here. 

1.5.2. The work of Tip. The recent papers of Tip U5JS5] is more closely related to 
this work. For the special case of the EM field in a linear absorptive dielectric, he 
has given a Hamiltonian formalism involving auxiliary fields similar to our "hidden 
string." His formalism made possible an analysis of energy conservation, scattering, 
and quantization [45] and led to a clarification of the issue of boundary conditions 
in piecewise constant dielectrics . Stallinga [32] has used this formalism to give 
formulas for the energy density and stress tensor in dielectric media. 

We do not rely on Tip's work, but the present work follows and parallels it to 
some extent. In particular, the present paper gives a general context in which some 
results of [45] [46] may be seen as special cases of results valid in a large class of 
linear dispersive Hamiltonian systems. 

1.5.3. Heat bath and coupling. The evolution equations (|1.251 fl~2l)]) describing the 
hidden string are identical to those of a so-called canonical heat bath as defined 
in [TH Section 2], [37] Section 2]. As the canonical heat bath as described in [TU 
Section 2], (37] Section 2] has naturally appeared in our construction, let us look 
at it in more detail. 

The Hamiltonian of our extended system (|1.18|) can be expressed as a sum of 
two contributions 

(1.52) R(U) = H sys (f7)+H str (t/), 
the system energy 

(1.53) H sys ((7) = ±\\Ku-T<pf H = i||/|| 2 , 

where / is the kinematical stress as defined by (|1.27[) . and the string energy 

2 , II Q / M|2 



(1.54) H str ([/) = i / H^sill;, + \\iK r -[s)r il 



ds. 



We conceive of H sys as the energy of an open system dynamically coupled to a "heat 
bath" with energy H str , described by the hidden string. 
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The physical concept of a heat bath originates in statstical mechanics. There gen- 
eral considerations show that for a system to behave according to thermodynamics 
it should be properly coupled to a heat bath. Dynamical models at the mathe- 
matical level of rigor were introduced, motivated and described rather recently, see 
[TBI Section 1], [T4J Section 2], |37j Section 2] and references therein. According 
to the references, based on statistical mechanical arguments, the heat bath must 
be governed by a self-adjoint operator with absolutely continuous spectrum and no 
gaps, i.e. the spectrum must be the entire real line K, and the spectrum must be of 
a uniform multiplicity. These requirements lead to a system equivalent to a system 
with the Hamiltonian H stv (ip, 9) as in (|1.54[> . t 14} Section 2]. 

Our construction of the unique extended Hamiltonian produces an auxiliary 
system with Hamiltonian in the form H str (^, 9) as in (|1.54|) and gives another way 
to obtain the canonical heat bath as a natural part of the conservative system 
extending a dissipative and dispersive one under the condition of its causality. In 
our Hamiltonian setting (| 1 . 52[) the coupling (Ku, Tip) can be classified as the dipole 
approximation, |37( Section 1,2], associated with a bilinear form. 

1.6. Organization of the paper. The main body of this paper has two parts. 
The first, comprising Sections [2] - |4j is essentially the physics part of the paper. It 
consists of a formal derivation of the quadratic Hamiltonian extension ($2]), con- 
taining all relevant physical details, followed by an application of the extension to 
TDD wave equations ($3j) with Maxwell's equations for the electro-magnetic field 
in a TDD medium considered as a detailed example (@. In particular, in $3] we 
write the extended Hamiltonian for a TDD wave system as the integral of a local 
energy density and derive expressions for the energy flux and stress tensor. We 
also derive general approximations for the time average of these quantities in the 
special case of an almost mono-chromatic wave. In 21 we specialize these formulas 
to the Maxwell equations. 

The second part, consisting solely of Section El is a more detailed mathemat- 
ical examination of the quadratic Hamiltonian extension. Here we give a precise 
formulation and proof of the main results leading to Theorem 11.11 with a rigorous 
analysis of the unbounded operators involved. 

The appendices contain supplementary material, including a.) a brief review of 
notation and spectral theory for operators on real Hilbert spaces, b.) a sketch 
of the proof of Thm. II. 2\ c.) a few examples illustrating the application of our 
construction to non-linear systems with linear friction and d.) a derivation of the 
symmetric stress tensor for a system with a Lagrangian density, used in Section [5] 



2. Formal construction of a Hamiltonian 

We begin by analyzing extended systems of the type outlined in (11. 18111. 22]) . with 
an unspecified symmetric operator valued coupling function ?(s). It is a simple 
matter to obtain, via a formal calculation, evolution equations of the form (|1.11|) 
for the reduced system. In this way, we get a relation between the susceptibility 
and the coupling ? — (|2.9|) below. As it turns out, the symmetry (|1.13| and power 
dissipation (|1.14|) conditions are necessary and sufficient for inverting (|2.9|) to write 
q as a function of \- 
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The extensions we consider are described by a vector coordinate U taking values 
in the extended phase space V = V © L 2 (M, H) © L 2 {R, H), with U eV denoted 

(2.1) U = I 9(s) 

\</>0), 

Recall that we interpret (j>(s) and 9(s) as the displacement and momentum density 
of an H- valued string, consistent with the equations of motion 24111. 27) , namely 

(2.2) d t u{t) = JK T f(t)+p{t) 

(2.3) d t 6(s,t) = #¥>(«,*)+?(*)/(*), 

(2.4) %>(*,*) - 6(s,t) 
with kinematical stress /, 

roc 

(2.5) f{t) = Ku(t)- / &aq{a)y{a,t) . 



Here we take <r to be an (as yet) unspecified operator valued distribution. 

Upon eliminating 9 from (l2.3U2.4p . we find that the string displacement ip follows 
a driven wave equation 

(2.6) {d?-d 2 s }vis,t)=,is)fit). 

Taking / as given, we solve (|2.6j) for ip with initial values limt^_oo (pit) = 0, 
limt-^oo dttp(t) = 0, corresponding to the string at rest in the distant past: 



(2.7) pis,t) 



1 ^.OO PS + T 

- j o dr J da,i*)fit-r) 



where we have tacitly assumed that f(t) is integrable. Recalling that / is related 
to tp by (|2.5|) . we obtain the following equation relating / and u 

(2.8) f(t) + - ds dr / da?( S )?(cr)/(i-r) = Ku{t) , 



2 . 

which is of the form of the generalized material relation (|l.llb[) with susceptibility 

i i-co rS+T 

(2.9) ^=o ds da,is),ia) . 

^ J — OO J S — T 

Thus, the reduced system described by u is a TDD system of the form (jl.llj) . 
with susceptibility given by (|2.9[> . To construct a quadratic Hamiltonian extension 
to (ll.ll|) it essentially suffices to write the string coupling ^ as a function of the 
susceptibility \ by inverting (|2.9p . Note that the r.h.s. of (|2.9[) is a symmetric 
operator, so the symmetry condition (| 1 . 1 3[1 is certainly necessary. As we will see 
the power dissipation condition (|1 . 14|) is also necessary, and together the two are 
sufficient. 

Note that (|2.9|) holds also for r < 0, with := — f , provided we replace x by 
its odd extension x°, defined in (|1.39|) . Differentiating with respect to r then gives 

(2.10) d TX °ir) = I [°° d«?(«){?(« + r)+?(«-r)} . 



If x(0+) 7^ then x° has a jump discontinuity at and (|2.10p . which holds in the 
sense of distributions, implies that <; includes a Dirac delta contribution at s = 0. 
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To understand the role of the PDC (11.141) here, let us suppose that 

(2.11) xM = limx( w + "?) 

exists and is continuous for to G R, as holds for instance if x G ^ 1 [0, oo). Then the 
PDC (fTTT4|) implies that 

(2.12) wlmxM > , 
which may be expressed as 

(2.13) wlmx(w) = -S^M, 
where 

(2.14) 5?(w) - / dre iwT x°(r) = 2i / dr sin(o;r)x(T) . 



To solve for we take the Fourier transform of (12.10p . which by (I2.13|) is 

(2.15) 2u;ImxM = \ + , 
with T(cj) = dse iws ?(s). Note that 

(2.16) ?(c)t = ?(_a;) , 

where •> denotes the Hermitian conjugate. Therefore (|2.15p is the same as 

(2.17) 2wIm X H = i^-w^-^t+fHfMt} . 

Clearly the r.h.s. is non-negative and we see, in particular, that ()2.9|1 implies the 
power dissipation condition (|1.14[) . (Once the inequality wlmx(w) > is known 
on the real axis, it extends to the entire upper half plane because Im£x(C) is a 
harmonic function. See (|5.32[ [5733)) below.) 

A solution to (|2.9p is not unique. However, there is a unique solution with <f(w) 
a non-negative real symmetric operator for each w, i.e., 

(2.18) ?H = = ?(w) T = ?(w) t = ?(-w), 
and 

(2.19) > . 

Indeed, under the symmetry condition (|2.18p . eq. (|2.17p simplifies to 

(2.20) 2wlmx(w) = ?(w) 2 . 

(This is consistent since Imx(w) is a real operator as we see from the formula 
Imx(w) = Jg 00 disin(wt)x(i)-) The unique non-negative solution to (|2.20p is given 
by the operator square root, 



(2.21) ?(w) = v /2wlmx(w) . 

Thus, a quadratic Hamiltonian extension of the system (jl.lip is given by (|1.18p 
ra'i/i i/ie coupling function <r given by Fourier inversion of the r.h.s. of (|2.2ip . i.e., 

(2.22) ?(s) = 77- / dw cos(wsV2wImxH , 
which is (|1.16p . 
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2.1. A minimal extension. The system with Hamiltonian (|1.18|) has a mechan- 
ical interpretation as strings coupled to the degrees of freedom of the underlying 
TDD system and provides a conceptual picture of the TDD dynamics in terms of 
absorption and emission of energy by those "hidden" strings. However, for calcula- 
tions and to describe the minimal extension, it is easier to work with a system in 
which the string displacement is replaced by its Fourier transform 

e iKS <p(s)ds. 

-OO 

To make the change of variables symplectic, we replace the momentum density 9 

by 

(2.24) 9{k) = — I e iKS 6»(s)ds. 

The resulting transformation of phase space 

/l 0\ 

(2.25) U h+ MU = ±T U, 

\0 T) 

with Tf{r£) — J R e 1KS f(s)ds, is a symplectic map — MJM T = J, since T~ x = 
■^J rT - Note that the Fourier transform T maps the symmetric space S(R,H), 
defined in (| 1 .48(1 . onto itself, since <p(s) — <fi(—s) implies that F4>(k) is real. Thus 
M. is well defined as a symplectic map on the reduced phase space V s (see f|1.47]> ). 
Correspondingly, we transform stress space by the orthogonal map 

AO ON 

(2.26) F ^ OF = -^T p. 

V° 



The map \T maps the anti-symmetric space A(M., H) (sec (|1.48|l ) onto itself, so O 
is well defined as an orthogonal map of the reduced stress space H s (see (|1.47p ). 
Together the two transformations amount to an symplectic/orthogonal isomorphism 
of the form (|1.49[) . and the impedance is transformed to 

(2.27) K ^ IC = OK,M-\ 

where 

(2.28) 
Here 

(2.29) Tip = I s(K)tp(K)dK, 





where fwas defined in (|2.2ip . 

The associated equations of motion, from the Fourier transform of (|2.2II2.5[) . are 

(2.30) d t u{t) = JK T f(t) 

(2.31) d t 9( K ,t) = -^ K 2 $( K ,t) + ^-s{K)f(t) 

Z7T ZTT 

(2.32) d t <p(K,t) = 2%e(K,t), 
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with 
(2.33) 



f(t) = Ku{t) 



1 

2t7 



<T(k)</?(k, t)dn. 



Combining (|2.31[) and (|2.32[) we obtain the Fourier Transform of (|2.6 
(2.34) d^(n,t) = -«V(M)+?(«0/(i), 

with solution 

rt sin(K(< - t')) 



(2.35) 



Clearly the resulting string displacement satisfies 

(2.36) <p{K,t) € ran (?(«;)) for every n E 



The same holds for the momentum density 9, since 

rt 



(2.37) 



0(M) = 8t<p(it,t) = 2its{k) 



cos(K(t - t'))f(t')dt' 



Thus, we may restrict the phase space to the Hilbert space 

(2.38) V = V ®S ? (R,H)®S ? {R,H), 
where 

(2.39) 5 f (R, H) = {/ S 5(R, H) : /(«;) e ran (?(k)) for every kg!}. 

We denote by Jo and /Co the restrictions of the symplectic operator J and impedance 
JC to Vo- Thus j7o still has the block matrix form (|1.22[) and /Co is defined by the 
r.h.s. of (|2.28[) for vectors U = (u,0,<j>) £ Vo- We consider the impedance /Co as a 
map from Vo to the restricted stress space 

(2.40) Ho = H ® S?(R, H) ® A ? (R, H) 
with (see (TOgl) ) 

(2.41) .4 f (R,i2") = {feA(J8L,H) : /(as) € ran (?(«)) for every At e M} . 



Clearly j7o,/Co give a quadratic Hamiltonian extension to (II. lip . We claim that 
the resulting extension to (| 1 .34[) is the unique minimal extension guaranteed by 
Theorem II .21 Indeed the generator has the expression 



/kjk 1 



(2.42) 



fcoJo^o — 



V 



/2tt 




-?=T 

V/27T 



At 







where [T T /](At) = ?"(«)/, by (|2.29p . One easily verifies there is no subspace of Ho 
invariant under /CoiTo/Cj and containing H © © (the range of Ih)- Thus: 

Theorem 2.1. There exists a quadratic Hamiltonian extension with (7i,/C 1 //C T ) 
the unique minimal extension of Theorem ] 1.2[ 

For the purpose of calculation it is sometimes useful to take the Fourier-Laplace 
transform (|1.15[) with respect to time, setting 

"(C) \ r oo / u(t) 



(2.43) 



6(k,0 



0(ac,*) I dt, ImC > 0. 
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(2.44) 


-Ku(0 


(2.45) 




(2.46) 




with 




(2.47) 


1(0 = 



(2-48) $(k,0 = - 5 -^ j ?(k)/(C) j 



(2.49) /(C) + — / ^— ^ K ) 2 dKf(() = Ku(Q. 



We obtain the system of equations 

= JK T f(C) 

= -i- K V(^C) + ^?M/(C) 

= 27r£(/t,0, 

1 f°° 

Ku ^~2ttJ ?( k )^K CMk. 
In particular (12.451 12.46|l together imply 

1 

\ 

which, with (|2.47[) . yields 

1 

2n «a - C 
This is suggestive of the identity 

< 2 - 50 > «o 'hL^e^ 2 ^ 

which holds for fas in (|2.2ip as we shall see in the proof of Theorem 15.41 below. In 
fact, (|2.50[) is a consequence of the Herglotz-Nevanlina representation for an (oper- 
ator valued) analytic function in the upper half plane with non-negative imaginary 
part (see [TJ Section 59] and [Ml Section 32.3]) or, what is essentially the same, the 
Kramers-Kronigs relations (see [HJ Sec. 62]). 

2.2. TDD Lagrangian systems. In many applications the phase space V decom- 
poses naturally as Vp®Vq (V p — V q ), with the symplectic operator in the canonical 
form0 

(2.5D j=(; - 1 

The two components of u = I ] <EV are "momentum" and "coordinate" respec- 



\1j 

tively. If the linear map K is block diagonal 



(2.52) K 



K p 
K q 



then the linear maps K p and K q can be thought of as follows: 

KpK p is the inverse mass (mobility) operator, and 

K^Kq is the stiffness (inverse flexibility) operator. 

Correspondingly, we suppose that K p is boundedly invertible, or at least invertible, 
as otherwise there are "infinitely massive" modes. The equations of motion are 

^Any symplectic operator can be written in the form 1 12. 511 by a suitable choice of basis for V, 
abut the subspaces Vp, q are not unique. See Lemma |A. 31 
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A system in this form has an equivalent Lagrangian formulation, with Lagrangian 

(2.54) Hq,d t q) = (j>,d t q)-h(p,q), 

where we express p as a function of dtq using the equation of motion for q, i.e., 

(2.55) p=[K^K p y 1 d t q. 

As K p is boundedly invertible, eq. (|2.55[) is unambiguous. Thus, the Lagrangian is 



(2.56) 



Hq,d t q) 



dtq 



h\K q q\\ 2 Ha 



where we have H = H p H q with K w g C(V W ,H W ), w = p, q. The trajectory 
q(t) may be obtained from the Lagrangian by noting that it is a stationary point 
for the action A([g(-)]; to, *i) = f t * L(q(t),d t q(t))dt. The Euler-Lagrange equation 
obtained by setting the variation of A equal to zero is 



(2.57) 



[K*K P ] l dU = -K^K q q, 



which is formally equivalent to (jl.lOal fl.lObp . with 



q 



For a Lagrangian system of this form, we generally make the physically natural 
assumption that an external driving force p(t) couples through the r.h.s. of (|2.57l) . 
That is the equation of motion is 

(2.58) [K^KvY 1 d 2 t q{t) = - K q q(t) + p(t) , 

with p(t) G V p = V q . This amounts to consider the time dependent Lagrangian 
L(q,d t q,) + {q{t),p{t)) y , i or Hamiltonian h(q,d t q) - (q(t), p(t)) V<l - 

For a TDD system (jl.llj) with Hamiltonian h of this form, the extended Hamil- 
tonian (081) . R(U) = ±\\ICU\\ 2 , is of the form 



(Kp -T p \ 



(2.59) 



K 



K q 




VO 



1 
d s J 



where T p , T q are the p and q components of the coupling operator T (see (|1.17|l ): 



(2.60) 



with momentum and coordinate string coupling functions ^ p and <j q respectively. 
Notice that the constitutive relation (12.511 turns into 



(2-61) / q = K q q - T q ip, f p = Kpp - T v ip, 

readily implying the following representation for the Hamiltonian 

H(C0 = \{\\kf H<l + \\U\\ 2 Hv } + \[ 
( 2 -62) = ±{\\K q q-T q ipf H ^ + \\K p p-T p <pf Hp } 



d s ip(s)\\ 2 H 



ds 



>( s )\\ 2 H + \\d sip (s)\\ 2 H 



where H = H a 
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We form a Lagrangian for the extended system, taking p and 9 as momenta, 

(2.63) L(q,<p,dtq,d t tp) = (p,d t q) + (6,d t <p) -K(p,q,6,ip). 
where we must write p and 9 as functions of dtq, dtf and (p, 

(2.64) 9 = d t ip, and p = [K%K P ] ~* d t q + K~%ip, 
using the equations of motion. Thus the above and (|2.63[) imply 



(2.65) L(q,<p,dtq,dt<p) = ~ [K^f d t q * + ( [K^f d t q,T p <p 

1 f°° 2 1 2 1 f°° 2 

+ 2j W d M s )\\ H ds - 2 W^-^Wh^ - 2 _/ ll^( s )llff ds - 
The Lagrangian form of the equations of motion, with driving force p, is 

(2.66) dt^K^K^dtqify+K- 1 j Sp(«)^( Sj «)d«} = -^ q T / q (<) + p(f), 

(2.67) flfo(M) = 9 S V(M) + Sq(«)/q(*) ~ C P («) [K- 1 ]' 1 d t q(t), 
with 

/•OO 

(2.68) / q (t) = K q q(t)- / c q ( S )^(s,t)ds. 



3. Local TDD Lagragians and conserved currents 

Many physical systems of interest are described by wave motion of vector valued 
fields, with the coordinate variable u a function of the position f £ K d (usu. d = 3) 
valued in a Hilbert space V. That is, phase space V = L 2 (R d ;V). Of particular 
interest are systems with a Hamiltonian that is an integral over R d of a density, 
whose value at a point r is a function of the field u{r) and its derivatives at the point 
r. In this section we focus on extended TDD Lagrangian systems of this type with 
u(f) — (p(r),q(f)) and the symplectic operator J in the canonical representation 
(123TD . 

That is, we take a system of the type considered in £12.21 and suppose the spaces 
Vq = V p and H Pyq are of the form 

(3.1) Vq = L 2 (R d , V Q ) and H w = L 2 (R d , H^), w = p,q, 

with Vq, JTp real Hilbert spaces. So the coordinate q is a vector field q(r) G Vb- 
We suppose further that the impedance operator K is of the form 

where repeated indices are summed from i = 1, . . . , d. For each r, -Kp(r) and if q (r), 
Yj(r), j = 1, ...,c? are bounded operators from Vq to H® and U q respectively and 
di = d/dp r This form covers classical linear elastic, acoustic and dielectric media. 
The system, in the absence of dissipation, is governed by a Lagrangian 

(3.3) Hq,d t q) = [ L(q(r),X7q(r),d t q(r))d d r, 
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with a Lagrangian density second order in the coordinate q(r) and first derivatives: 



(3.4) L(q(r),Vq(r),d t q(r)) 
1 
2 



K p (r) T ] d t q{f) ~ \\K q (r)q(r) + Y(r) • Vg(r)||' 



By a suitable choice of K Pyq and Y we can obtain any Lagrangian density of the 
form L = T(dtq) — V(q, Vg) with T and V homogeneous of degree two. 

Now suppose there is time dispersion and dissipation in the system so that the 
equations of motion and material relations according to (II. lip and (|3.2[) are 



(3.5) 



0, 



p(r, t) 
q(r,t) 



-if q (r) T / q (r, t) + ft { Y(r) T / q (f, t)} 



(3.6) 



f P (r,t) 



y(T . i *(fp(f>t-T) 



dr = 



K p p(r,t) 



^K q q(r,t)+Y l (r)d l q(r,t) / 

with x( T i r) a B(Hp © H®) valued susceptibility functionQ The string coupling 
operators constructed above then fiber over R d in the same way 

(3.7) [<r p (s)p](r) = s p (s,r)p(r), [s q {s)q] (r) = <r q (s, r)q(r), 

and the extended Lagrangian (|2.65p is the integral of a Lagrangian density: 



(3.8) L(Q, d t Q) = 

where Q(r) — (q(r), ip(s,r)) and 

(3.9) L(Q(f),VQ(f),ftQ(r)) 
1 



L(Q(f),VQ(r),ftQ(r))d d r, 



[K^(f)] 1 d t q(f)f H0 + ^ [X p T (r)] 1 d t q(f), J™ , p (s, r>( S , f)ds^ 



H" 



\dMs,r)\\ HOfBHO ds 



dM^OIIffOg^O d.s 



K q (r)q(r) + Y(r) ■ Vq{r) - / ? q (s, f)(p(s, r)ds 



HI' 



Remark: ip is an element of L 2 (R, H) with H = H p © H q = L 2 (R d , H° © H°). We 
identify L 2 (R, H) with L 2 (RxR d , ff°©ir°), writing cp as (s, f) i-> (p(s, r) € H°®H°. 

In Appendix[D]we recall some basic constructions for a system with a Lagrangian 
density. In particular, we obtain suitable expressions for the energy flux vector and 
the stress tensor of a homogeneous system. In this section, we apply the expressions 
derived there to the extended TDD Lagrangian (|3.9[) 



This form for the susceptibility precludes spatial dispersion, which would involve integration 
over r on the l.h.s. of l|3.6| l. Our general construction works in the presence of spatial dispersion, 
but the extended Lagrangian is non-local making it difficult to give a meaningful definition of the 
energy density and stress tensor. 



22 



ALEX FIGOTIN AND JEFFREY H. SCHENKER 



3.1. Energy density and flux. Because L does not depend explicitly on time, 
the total energy, which can be expressed as the integral of a density (see (|D.4j) ) 

(3.10) £ = [ H(r,t)d d r, 

is conserved (in the absence of an external driving force). The value of the total 
energy £ is, of course, just the extended Hamiltonian H evaluated "on-shell," at a 
field configuration evolving according to the equations of motion. 
We can express the energy density H as a sum of two contributions 

(3.11) H(r,i) = H sys (r,t) + H str (r,i), 

which we interpret as the energy density of the TDD system and the heat bath, as 
described by the hidden strings, respectively. Here 

1 f°° 1 f°° 

(3.12) H str (r,t) = - / \\d s ip{s,r,t)\\ HO(BHO ds + - / \\9(s,r,t)\\ H0(BH0 ds, 

L J -oo p q Z J_ QO p i 

with 6{s,r,t) = d t tp{s,r,t), and 

(3.13) H sys (r,t) = \ \\UrM 2 H o + \\\U^t)\\ 2 H », 

with /q jP the coordinate and momentum parts of the kinematical stress, 



(3.14) / q (f) = K q (r)q(r) + Y(r) • Vq{f) - / ? q (s, f)ip(s, r)ds, 

J — oo 

and 

(3.15) f p (r) = [ifj(r)] d t q(r,t) = K p (r)p(f,t)- / ?p (a, r>(s, f, t)di 



as follows from (|2.5p . (|3.2|) and (|3.6p . By (|3.1lH3~TS|) . H sys includes the interaction 
energy between the system and the hidden strings. 

The energy density, expressed in canonical coordinates, H(r, t) = H(g(r, t), 
Vq(f,t), Lp(-,r,t), p(f,t), 0(-,r,t)), is also the Hamiltonian density. The equations 
of motion can be recovered by variation 

(3.16) d t q{r, t) = — (r, t) = K p T (r)/ p (r, t), 

(3.17) d t p(f,t) = +9^} (r,t) = -K^(r)f q (r,t) + diYj^f q (r,t), 

(3.18) d t 6(s,f,t) = --—(r,t) 

Sip(s) 

= d 2 Ms, f, t) + s)f p (f, t) + <^(r, s)/ q (r, t) 
8Y\ 

(3.19) d t y{s,r,t) = j^(r,t) = 6(s,r,t). 

oo(s) 

If the system is driven by an external force p(r, t) £ Vq, we replace (13. 17)) by 

(3.20) d tP (r,t) = -K^^f^^ + diYJirJf^+p^t). 

When the driving force vanishes, the total energy is conserved and the energy 
density satisfies a local conservation law 

(3.21) d t H(r,t) + diSii^t) = 0, 
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with S the energy flux vector, an expression for which is derived Appendix |D] For 

the present case the energy flux is (see (|D.5jl ) 

(3.22) 

S,(r,t) = -(d t q(r,t),YT(^f q (r,t)) Vo = - (K? (r) f p (r, t), Yj(f)f q (f, t)) Vo . 
With a non-zero driving force, one has 

(3.23) d t H(f,t) + diSi(f,t) = {d t q(r,t),p{r,t)) Vo 
in place of (|3~2Tj) (see Theorem EH • Thus 

(3.24) d t £ = [ (d t q(r,t),p(f,t)) Vo d d r, 

consistent with our interpretation of (dtq(f, t), p(f, i))y as the power density of the 
external force. 

We conceive of the four vector fields p(f, t), q(r, t), f p (f, t), and / q (r, t) as spec- 
ifying the "state" of the reduced TDD system. By (|3.22j) , the energy flux at time t 
is a function of these fields evaluated at time The energy density depends in a 
more essential way on the configuration of the hidden strings. Nonetheless, we may 
use (|3.23[) to give a definition of the energy density intrinsic to the TDD system by 
writing it as integral over the history of the system, namely 

(3.25) H(r,*) = f {-5 i S i (r,t , ) + (5 t 9(r,t'),p(f;t , ))vo}cit / ! 

J — oo 

if the energy density was zero at t = — oo, i.e. the system and medium were at rest. 

3.2. Homogeneity, isotropy, wave momentum and the stress tensor. Sup- 
pose the extended TDD system has a Lagrangian density (|3 - 9(1 which is homoge- 
neous — invariant under spatial translations. Then the total wave momentum P 
is conserved, and there is a corresponding conserved current, the wave momentum 
density p(r, t), analyzable using Noether's Theorem (see [13 Section 5.5] )Q If the 
Lagrangian density is also isotropic — invariant under spatial rotations — , then the 
anti-symmetric tensor of angular momentum about the origin Mi j is conserved!! In 
appendix [D] following [3], we recall the formulation of the symmetric stress tensor 
T and wave momentum density p for a homogeneous and isotropic system. 

To say that the Lagrangian in (|3.9[) is isotropic, we must specify how q and 
ip transform under rotations. Thus we suppose given representations T and T w , 
w = p, q, of the rotation group SO(d) by orthogonal operators in B(Vq) and B(H^), 
w = p, q, respectively, so that under a global rotation of the coordinate system 
about the origin, 

(3.26) r' = R-f 

with R £ SO(d) an orthogonal matrix, the fields q and tp transform as 

(3.27) q'{r') = T(R)q(R- 1 f') 

^This is a consequence of the absence of spatial dispersion. By adding terms involving \/<p to 
the Lagrangian, we could extend the above set up to systems with spatial dispersion, resulting in 
an energy flux with non-trivial contributions from ip. 

^We follow |30] in using the term "wave momentum" for the conserved quantity associated to 
translation invariance. This avoids confusion with "canonical momenta," the variables p and 8. 

dimension d = 3, the usual angular momentum pseudo-vector t is obtain from Mi j as 
&i = ei,j,fcMj,fc with €i j k the fully anti-symmetric symbol with £1,2,3 = !• 
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and 

(3.28) <<f(a,?') = ( T *^ T ^ R) ^( s ,R-^r'). 

(Recall that if is H p © Hq-valued.) The representations T and T Pi q are specified by 
families of skew-adjoint operators g^j and gf^, w = p,q, on Vq and H^, w = p, q, 
respectively, i, j — 1, ...,d. A rotation R = e u , with uj an anti-symmetric matrix, 
has representatives (see appendix |D|) : 

(3.29) T(e w ) = e^' 9 ", T w (e") = e* Ut - s9 *-j (summation convention) . 
In addition to being skew-adjoint, the operators gij and g?'^ satisfy 

(3.30) gl = -g]., 

and AA',3' - al'j'dtj = -ki'9j, r + hf9j,i> + hi'9h> - hraii" 

Since Vq and H p q are finite dimensional in our application to the Maxwell equations 
below, we assume gij and gfj are bounded for simplicity. 

Definition 3.1. The Lagrangian density (|3.9[) is homogeneous if K p q (r), Y(r), 
and ? p ,q(r, s) are independent of r, and is isotropic if 

(3.31) K w g lyj = g™ 3 K w , w = p,q 

(3.32) Y; ,/,.., = g^Y, ■ ■i.j.Y, ,),j,Y, 

and 

(3.33) w(s) (^^ ff qJ = S*jW(s), w = p,q. 

Remarks: i.) In Appendix [Dl we give more general definitions ID . 1 1 and ID . 2i which 
are consistent with 13.1 1 ii.) The last two terms on the r.h.s. of (|3.32|) result from 
the fact that Yj appear coupled with a spatial derivative in the Lagrangian (|3.9[) . 
iii.) Recall that ? w (s, r) : H p © _ff q — > _ff w . 

Theorem ID. II below gives the following expressions for the wave momentum 
density p and stress tensor T, expressed here in canonical coordinates: 

(3.34) = (d iq (r,t),p(r,t)) Va 

+ / (di¥>(s,r,t),6(s,r,t)) H o mH<J ds - dj$ it j(r,t), 



(3.35) T 4J (r,i) 

= - (d t q(r, t), Yj/ q (r, i))^ - ^-L^, i) + 5 t $ tJ + d^,,- k (f, t), 
where / p . q are defined in (|3.14l |3~T5|) . L is the Lagrangian density 

(3.36) L(r,i) = <KT/ p (r,t),p(r,t)) H0 -i||/ p (f,i)||^-i||/ q (r,t)||^ 

+ ~ / \\0{r,s,t)\\ 2 H „ds- - / ||a s v3(s,r)||^o e// o ds, 
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and 

(3.37) $ id (?,t) = ^(g hJ q(r,t),p(r,t)) v0 

9 l 3 i)<p{s,f,t),e(s,r,t)) ds, 

(3.38) ttu.fcCf.t) = -l(g^q(r,t),Yjf q (r,t)) Vo 

+ \ (gj,kQ(f, *), Yjf q (r, *)>y + \ (9iM^ *), Yj/ q (r, t)) % . 
As the system is homogeneous, the total wave momentum 

(3.39) P = / p(r,t)d d r 



is conserved in the absence of a driving force, and the wave momentum density 
satisfies the local conservation law 

(3.40) d tPl (r,t) + d J J tJ (f,t) = 0. 
With a driving force, (|3.40j) is modified to (see Theorem ID. 2p 

(3.41) ftPi(»*,*)+ajT w (f,t) - {d iq (f,t),p{r,t)) Vo , 
and thus 

(3.42) dtP = / (^(r,t),p(f,0)y o d< i f. 

Due to the term dt&ij in the definition of the stress tensor (|3.35|) . the driving force 
also modifies T. If T? -(r, t) is the stress tensor with p(f, t) = 0, then 

(3.43) Tij (r, t) = (r, t) + i ( Siii g(f, i), p(r, t)) Vo . 

Proposition 3.1. If the Lagrangian is homogeneous and isotropic then the stress 
tensor Q3.35P can be written as 

(3.44) T id (f,t) = 5 {Wy + W w (fit)} + i (sjjg^t),^,*))^ - ifjL^t), 
where p(f, t) £ Vo is the external force and 

(3.45) Wij(f,t) = -(d iq (r,t),Yjf q (r,t)) Vo +d k (g hk q(r,t),Yjf q (f,t)) Vg , 

with f„, f p and g id as defined in (|3.14p . (|3.15[) . and (|3.271f3~2"9")) . In particular, the 
stress tensor is symmetric if p(r, t) — 0, and depends on the state of the strings 
ip(s,r,t) through the variables f q and f p and the Lagrangian density L. 

Proof. By the remark preceding the theorem, it suffices to consider p(r, t) = 0. 
Using the equations of motion ()3 . 1 61 T5TT5|) . the definition of isotropv l3.11 and the 
skew-adjointness of the rotation generators g\ j , one may calculate that 

(3-46) = l(9ijq,d k YJU) v .+~(g? tj Y k d k q,f <l ) Yo 

= d k^(g l , 3 q,Ylf q ) Va - ({Y k g itj - gfjYk} d k q, f q ) H ^ 

= d k \( 9l ^Ylf q ) Va Ud jq ,Yjf q ) Vo + Ud iq ,Yjf q ) Vo . 
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Combining this expression with the definition (|3.35|) of T gives (|3.44[) . □ 

Like the energy fiux S, the tensor field W at time t is a function of the fields 
q(r,t), p(r,t), f q (r,t), and f p (r,t) specifying the state of the reduced TDD sys- 
tem. In particular, the off-diagonal terms of the stress tensor depend only on the 
instantaneous state of the reduced TDD system. The diagonal terms depend on 
the Lagrangian density L, requiring a more detailed knowledge of the state of the 
hidden strings. However, as with the energy density H, we may write L in terms of 
the history of the underlying TDD system. To this end, we rewrite L as 

(3.47) 

rOQ 

L(r,i) = {p{r,t),d t q{r,t)) Va + f (0(s, r, t), d t <p(s, r, t)) H0@H0 ds - H(r, t) 



= (K p p{r,t),f p {r,t)) H o + \\6(s,r,t)\\„ 0(BH0 ds-H{r,t), 

J — oo p q 

using the equations of motion (I3.16ll3~19|) . Based on (|3.11hl3~T3|) we introduce 

(3.48) L(r,t) = L sys (f,i) + L str (f,i), 
with 

(3.49) L sys (r,i) = (p(r,t),d t q(r,t)) Vo -H sys (r,t) 

= (K p p(f, t), / p (f, t)) H o - \ ||/ p (f, t)f H0 - \ ||/ q (r, t)\\ 2 Ha , 

poo 

(3.50) L str (r,i)=/ (0(s, r, t), d t <p(s, r, i)) ff o eiJ -o ds - H str (r, t) = 



1 f 1 

- \\&(s, r, t)\\ H o mH o ds- - I \\d 8 (p(s,f,t)\\ d H a Bm ds. 



We use the solution (|2.7[) to express 6(s, r, t) as 
(3.51) 6(s,r,t) = d t <p(s,r,t) 

1 

2 



t 

dt' { c (s + t-t')+ ? ( fl - t + t')} f(r, 0, 



where ? = (? p ? q ) , and / = ( j P ) • Then by the definition (|2T2"2"j) of see (f^TTU]) . 



OC 

2 



(3.52) / \\6(s,r,t)r m@H ods 

(f{f,h), [d T X°] {tl - t 2 )f(f f ,t 2 )) H o mHO dhdt 2 



— oo 

t 



1 



oo </ — oo 
<t r t 



(f(r,h), [d T x°] (2t - ti - t 2 )/(f ? ,t 2 )) ff o ffiff0 dt!dt 2 



oo — oo 



where x° is the odd extension (|1.39p of the susceptibility. Using (|3.25p to express 
H and (|3.52[) to express the corresponding term in (|3.47|) . we obtain an intrinsic 
definition of the Lagrangian density, and hence the stress tensor T, as function of 
the history of a TDD Hamiltonian system. Writing 

(3.53) p(r,i) = f {-9 j T iJ (f,0 + Mf,0>p(^i'))}dt , ) 

J — OO 
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we obtain a similar expression for the wave momentum density. 

3.3. Brillouin-type formulas for time averages. As we have seen, to express 
the energy density and stress tensor of the extended system in terms of the fields 
P) 9) /pj /q we must introduce integrals over the history, like (|3~25|) and ([3752]) . How- 
ever, it is often useful to have an approximate formula involving the instantaneous 
state of the TDD system. A well known example is the Brillouin formula for time 
averaged energy density stored in a dielectric medium (see [22j §80] and JO] below). 

Taking inspiration from the Brillouin formula, we consider here an evolution of 
the underlying TDD system which is approximately periodic with frequency uj/2tt. 
That is, we suppose that 

(3.54) g(r,t) = Re {e^ g (r, t)} , 

with g = p, q, f p , / q , or p. The various functions w = p , q , / p;0 , / q; o, Po are 
supposed to vary extremely slowly over time scales of duration 1/w, and may take 
values in the complex Hilbert spaces CVq, CH® . This evolution describes a carrier 
wave of frequency lo/2tt, which is slowly modulated in phase and amplitude. 

To quantify the notion that the functions go vary extremely slowly on time scale 
l/u>, we assume the Fourier Laplace transforms, 

/•OO 

(3.55) g {r,0 = / e^g (r,t)dt, Im( > 0, 



for g = po, q , / p;0 , / q;0 , or p , satisfy 

(3.56) \\go(r,()\\ < const. Wo V(ICIM), ImC > 0, 

with ip a fixed rapidly decaying function and ujq << oj. Thus 8 — loqIlo is a 
dimensionless small parameter which measures the slowness of the functions go- We 
are interested in asymptotic expressions for various quantities as 6 — > carried out 
to order 5 and shall neglect contributions of size o(S). Throughout the discussion 
the carrier wave frequency ui is fixed, so 5 oc ujq. (Recall that o(S) denotes any term 
with o(S)/S — > as S — > and 0(6) denotes a term bounded by const, x 5.) 

We use the notation a w b to indicate that a — b = o(8) and say that a is negligible 
if a 1=3 0, i.e., a = o(6). For example dfgo(r, t) is negligible for each go, since 



i 

< const. / v 2 ^ ( | v\ /ujq )dv/ujQ — O (o;q ) 



(3.57) \d?g (r,t)\ 



by ^M- Similarly (d t g (r,t)) 2 » 0, d?g (r,t) » 0, etc. 
We also write the string fields in the form (|3.54|1 . i.e., 

(3.58) 8(s,r,t) = Re {e~ [uJt 9o{s, f, t)} and ip{s,f,t) = Re {e~ iuJt ip {s, r, t)} 

However, it is convenient to use the formulation of I 



(3.59) fr r ?ri = r e iKS ( M iS 'flf) ds, 



<p(n,r,t)J V V>{s,f,t) 
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involving the Fourier transform of the string variable s. Then (13.58| implies 

(3.60) 6(K,r,t) = ^{e- iu %(K,r,t) + e iu) %(-K,f,t)*], 

and <p(K,r,t) = ^ {e^*^*, r, t) + e ibJt ^{~n, f, <)*} , 

where •* denotes complex conjugation. 

The string equations of motion (|3.18[ 13 . 19[) imply the following for 8q and ipo : 

k 2 ~ 1 

(3.61) d t o (K,f,t) -iuj0 o (K,r,t) = - — (p (K,r,t) + —s{n)fo(r,t), 

Ztt Zir 



(3.62) 



d t (p (K,r,t) -iutp (K,r,t) = 2tt9 q (k, r, t) 



The solution to (|3~6T1 13~62|) with O 



with ?(k) = (%(k) ? q («0) and f Q = (j Pl °^j ■ 
and ipo vanishing as t — * — oo is expressed, by the Fourier inversion formula 

(3.63) 9 (K,r,t) = ?(r,K) ' 



(3.64) 0o(re,r,t) = ?(r, k) 



4tt 2 
1 

2^ 



et _ il/4 -i(w + z/ + ie) 7 
e et wt — — ^ r-r^fa{r, v + \e)dv 



n et—iut 



k 2 ~ (lo + i/ + ie) 2 ' 
1 



, k 2 — (w + i/ + ie) 2 ' 

with e > arbitrary. 

The string energy density H str (r, i), as given by (|3.12|) , may be written 

,_2 



(3.65) H str (r,i) = i 



27r||e( K) r,t)|| 2 + — ||^(«,f,t)|| 2 

Z7T 



dn. 



Due to the dissipative dynamics of the reduced system, we expect a steady accu- 
mulation of energy in the string degrees of freedom. That is, H str (? ; ', t) should grow 
steadily until the work done by the external force p is completely dissipated. Thus 
H str should depend quite strongly on the history of the system. Thus we consider 
the rate of dissipation of energy to the strings, the power density 9tH str (r, t). 

On time scales of order 1/w, the power density <9tH s t r (r, t) may fluctuate wildly. 
To eliminate these fluctuations, we consider the time averaged power density 



(3.66) 



d t H str (r,t) = - f p(T/<i)d t H stI (r,t-T)dT, 



where /? is a fixed Schwarz class function with J /3(r)dr = 1 and a is a time scale 
much larger than l/u> but sufficiently short that /o varies slowly over intervals of 
length a, i..e 1/w < c < 1/wo. To provide for that with fixed carrier frequency ui 
and 6 = ^ -> we take 



(3.67) 

with < e < 1/2, readily implying 



1 



a cx S e x uj q 6 



(3.68) 

(Recall that 1 /luq is the time scale for /o variation and we consider the limit cjq — > 0.) 
We also assume that 

POO 

(3.69) / Tf3(r)dT = 0, 
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as holds, for instance, if (3 is symmetric about zero. Then given a slowly varying 
quantity Q{t), for which 

(3.70) Q(t - r) = Q(t) - rd t Q(t) + t 2 0(5 2 ), 
we have 

(3.71) Q(t) = Q(t)-d t Q(t)~ r/3(r/c7)dr + - / 0(^ 2 )r 2 /3(r/a)dr 

= Q(t) + 0(5V) « Q(t), 
since S 2 cr 2 = S^^w' 1 = o(S) for e < 1/2. 



Proposition 3.2. TTie time averaged power density of dissipation, <9 t H str , /ias i/ie 
following expression, to order o(S), 

1 ' 



(3.72) d t H stI (r,t) « (/„(r,t), W Imx(r»/ (r,t)) Cfl 



+ Im (d t f (r, t), d u wlmx{r, w) f Q (?, t)) CH 

1 

2 



+ o^t (fo(r,t),d^Rex(f,oj)f (r,t)) CH 



Remarks: 1.) The inner product (»,»)cjj denotes the complex inner product in 
Cif, linear in the second term and conjugate linear in the first. 2.) In general, the 
last two terms on the r.h.s. of (|3.72j) are 0(5). However, the first term is of order 
0(1) and is non-negative by the power dissipation condition. To first order, energy 
is dissipated at a steady rate governed by the size of the wlm£(r, to): 



(3.73) d t H stx (r,t) = -(f (r,t)^lmx(r,uj)fo(r,t)) CH + 0(S). 

Proof of Prov. HOI By Q3.65P and (|3.66p . the time averaged power density is the 
sum of two terms, which may be approximated as follows 

/>00 P oo 

(3.74) -/ /3(r/a) d t \\e( K , r, t - t)\\ 2 &k&t w - / i)i\\tt n {h\ v. /)||'"M/v. 



(3.75) -!- / /3(t/(t) / K 2 9 t ||^(K,f,i-T)|| 2 dKdr 



47TCT 



i / « 2 9t 1 1 («, t% *) 1 1 2 d«. 



On the r.h.s. 's of Q3.74[ I3.75|) we have dropped terms with the rapidly oscillating 
factor e ±2lllJ ' , as their time average is smaller than any power of 6 as can be seen by 
repeated integration by parts. Furthermore we have dropped time averaging from 
the remaining terms, by (|3.71|) . since we will show that J ||0o|| 2 dK and / K 2 ||^o|| 2 dK 
are slowly varying in the sense of (|3.70j) . 
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Let us first sketch the integration by parts argument allowing to neglect the 
terms dropped. We focus on a single term missing from the r.h.s. of (13.74| . namely 

poo pOO ^ 

(3.76) — / 0(T/a)d t e 2i "(*- T )(0oKr,i-r)*,0 o (-«,r,t-r)) d K dr 
4c 7_™ 7^ \ / ch 



— OC 

oo 



= j— 2 I e^Wpir/a) (e (K,f,t-rr,d (-K,r,t-T)) d^dr, 

where we have integrated by parts once. Although we have gained a factor of 
I/a, this does not yet show this term is small, because ||0o|| could be as large as 
1/5 oc cr 1 / 6 due to the large amount of energy absorbed by the strings up to time 
t. However, using exp(2iuit) = (2iw) _ ™i9™ exp(2iw£), we may integrate by parts as 
many times as we like. Thus for any n, the r.h.s. of (|3.76p equals 
(3.77) 

/OC f p OC ^ 

^ e^-r) d n ^ {t/u) j ^ ( 6o{Ki fi t _ T y,6o(-K, r, t - t)) ch d.j dr. 

Each t derivative acts either on (3 1 or on (8q, #o)- In the first case, we gain a factor 
of l/cr = lo5 £ and in the second case a factor of S. Thus this term is 0(5 ne ) and, 
as n is arbitrary, smaller than any power of 5. The other terms missing from the 
r.h.s. 's of (|3.74[ [3775)) — there are three in total — may be estimated similarly. 

To approximate the two terms on the r.h.s. 's of (|3.741 13.T5[) . we use the repre- 
sentations p. 631 [3TM|) for #o and ipo- For instance by (|3.63|) we have 

/oo 
d t \\e ( K ,f,t)\\ 2 d K 
-oo 

g ei{vi _„ 2 _ 2ie)t {u_+vx -ie){<jj + v 2 + ie) 



32^ 3 JJJ ' ( K 2 -(o; + i/i-ie) 2 )(K 2 -(w + i/ 2 + ie) 2 ) 



x (fo(r, v\ + ie), s{r, K) 2 f (r, i>2 + ie)) dz/idz/ 2 d«;. 

\ I CH 



CH 

Interchanging integrals to perform the n integration first, we compute 

f°° 1 

(3-79) / — — ■ , — ■ — ?(r, nfdn 

J-oo ( K - (w + z/i - ie) 2 ) (k- - + v 2 + ley) 

1 



(w + v i — ie) 2 — (to + V2 + ie) 2 
1 



• 2 — (a; + v\ — ie) 2 k 2 — (u + i/ 2 + ie) 2 



<f(r, k) d«; 



(2w + vi+ V2){v\ - V2 — 2ie) 
by (|2~50|) . Therefore, taking e -> 



(3.80) J / d t \\d a ( K ,f,t)WM K 



2n 

T^{x(r,uJ + vi +ie)* -x(r,w + z/ 2 + ie)} . 



2 



i6tt 2 yy 2w + 2/1 + 1/ 2 

x (/o(r>i), {x(r,w + i/i)* - x(r, w + i/ 2 )} /o(r> 2 )) AviAv 2 . 

\ I CH 
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Expanding to first order around V\ = v 2 = we have 

(3.81) i ( " + " lK " + V2) {X(r, » + vxT - X(r, u, + v.)} 

ZUJ + V\ + v% 

Thus 

(3.82) - / 9J0 o ( K ,r,t)|| 2 d K 

~ 7S (/o(r,t),wImxtr,u;)/o(f,t)) CH + Im (c't/o(^t),9 w wImx(r,w)/o(r,t))cH 



+ o^t (/o(r, f), w5^Rex(r, w)/ (r, t)) CH 



In a similar way, by (|3.64[) . 

1 f°° 

( 3 - 83 ) ^ _/ « 2 9 t ||^o(K,r,t)|| 2 d K 

3 te + i ( 1/ i- ly 2-2ie)t_ 



32tt 3 JJJ {k 2 - (w + vi - ie) 2 ) (k 2 -{uj + u 2 + it) 2 ) 

x (M^, *>i + ie), ^X^, K) 2 fo(r, v 2 + ie) ) d^id^drt. 

The approximation 
1 f°° 

(3.84) / K 2 cy<A,(K,r,t)|| 2 dAC 

~ I j (/o(^i),wImx(r,w)/o(r,t))c# + Im {d t fo(r, t), d u ulmx(r, w)/ (f, t)) CH 
+ \d t (f (f,t), -d a u 2 Rex(r,u})f (?,t) 

Z \ w / CH _ 

follows, since 

/•«> K 2 

(3.85) / — -5 7 . . „ 7 — rrr ?Tr, «) 2 dK 

7-oo (k 2 - (w + ^1 - ie) 2 ) (k 2 - (w + v 2 + ie) 2 ) 

1 



(cj + v\ — ie) 2 - (uj + v 2 + ie) 2 

(w + vi- ie) 2 (w + ^2 + ie) 2 



k 2 — (u> + v\ — ie) 2 k 2 — (uj + ^2 + ie) 2 
2tt 



<f(r, «;) 2 dK 



(2a; + v x + ^2X^1 — 1>2 — 2ie) 
x {(a; + i^i — ie) 2 x(^, w + u\ + ie)* — (cj + i/ 2 + ie) 2 x(r, uj + v 2 + ie)} , 

again by (|230| . 

Combining (|3~74l l3~75l EOWl l3~84l) we obtain (j3~72l) . □ 
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If there is no dissipation at r at frequency to, so 

(3.86) Imx(r, lo) = (zero dissipation at lo), 

then the string at r does not effectively absorb energy at frequency lo. We expect 
the total dissipated energy to fluctuate but not grow, so there should be a formula 
similar to (|3.72[) for the time average of H str . In this case, only the third term 
contributes to the r.h.s. of (|3.72|) . This term is a total derivative, suggesting the 
approximation 

(3.87) TU(f, t) = \ (/o(f, t), d^Rex(r, Lo)f (f, t)) CH + 0{8) (Im£(f, u) = 0). 

Using the methods of the proof of Prop. 13.21 one may verify that (|3.87[) is indeed 
correct. By similar arguments, we find that the time average of the system energy 
H sys (r, t), defined by (|3.13p . satisfies 

1 

4 



(3-88) H sys (f,t) ^-(f (f,t)J (f,t)) CH . 



whether or not there is dissipation at lo. Combining (|3.87| and (|3 . 88[1 we obtain: 



(3.89) H(r,t) = 7 {f (?,t),d u [uj(l + Rsx(f,u;))]fo(r,t)) CH + 0{S). 



Proposition 3.3. If there is no dissipation at frequency lo at r then the time 
average of the energy density at f satisfies 

1 

4 

For Maxwell's equations in a TDD dielectric, the above formula (|3.89|) reduces to 
the classical Brillouin formula for the energy density, [HJ §80]. (See (|4.71[) below.) 
Thus, (|3.72p may be viewed as an extension of the Brillouin formula to frequencies 
with dissipation and to arbitrary TDD Hamiltonian systems. 

Similarly, we may consider the time averaged string Lagrangian density 

I r°° If 00 

(3.90) Ur(r,t) « - / \\0o(s,r,t)fds-- ||c^o(s,r,t)|| 2 ds, 

where L str denotes the quantity 

(3.91) L str (f,i) = - / ||<?( S ,r,i)|| 2 d S -- / \\d s ip(s,r,t)\\ 2 ds. 

* J — oc " J-oo 

First note that by combining (|3.821 13.84| we have 

(3.92) d^(r,t) « ~d t (f Q (r,t),Re X (r,uj)f Q (r,t)} CH . 

This approximation holds (to order o(5)) whether or not the dissipation wlmx(a;) 
vanishes at frequency lo. Both sides are total time derivatives and, in fact, we have: 

Proposition 3.4. The time averaged string Lagrangian density satisfies 

(3.93) UAf,t) = -^(Ur,t),Re X (r,Lo)fo(r,t)) CH + O(S), 

and the total Lagrangian density L satisfies 
(3.94) 

L(r, t) = ^Re (/ p;0 , K pPo (r, t)) CH - i (/ (r, t), (1 + Rex(r, w)) f (r, t)) CH +0(6). 
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Proof. Eq. ([3J4]l follows from ([3^3]) using the expressions (j3~48l [3~50]) for L, the 
approximation (|3.88|) for the system energy density, and further integration by parts 
to obtain the first term on the r.h.s. from the corresponding term in (|3.49p . 
To verify (|3.93[) . we follow the proof of Prop. [3~2l to obtain 



(3.95) L str (f,t) 
1 



e +i(i/i-w 3 )t. 



167r 2 J J 2lo + v\ + v 2 v\ — v>2 — iO 

x {(fo(r,vi), [{uj + v-i)(u + v 2 ) - (lu + vx) 2 ] x(r, u + vi)*fo(r, v 2 )j 



(/o(r, vi), [(w + vx)(u + v 2 ) + v 2 ) 2 ] x(r,u> + v 2 )fo(r, ^a)) ^|d^id^ 2 



CH- 

1 



167T 2 J J 2ui + V\ + V 2 

x (fo(r,vi),[(uJ + i'i)x(r,u} + i'i)* - (u + v 2 )x(r, u + V2)} fo(r, V2)) d^idz/ 2 - 

\ / CH 

The key point is the cancelation of l/(Vi — v 2 — iO) by V\ — v 2 resulting from 
+ v\){u> + v 2 ) - (w + f2) 2 = -(w + fi)(fi - fa), 

(3.96) 

(a; + + 2/2) — (w + fi) = (w + ^2)^1 - v 2 ). 
Up to terms of order 0(S) we may replace 



1 

(3.97) - — ■ ■ (w + z/i)£(r ? ,o; + ^1)* - (w + v 2 )x(r,uj + v 2 )} 

Luj + v\ + v 2 

by its value at v\ = v 2 = 0, which is Rex(r, ui) resulting in (|3.93[) . □ 

Our main interest in (|3.94|) is in approximating the time averaged stress tensor. 

Proposition 3.5. For a homogeneous and isotropic system, the time averaged 
stress tensor satisfies 



hie {W, ,:„( /••• t) + W .„,,,(,••. t)} + h 

-S itj L(r,t)+0(S), 



(3.98) J lJ (f,t) = -Re{W lJ;0 (f,i)+W J - l; o(f,i)} + -Re( ffljgo (r,t),Po(r,t)) CH 



where L(r, t) is given by (|3.94|) and 
(3.99) 

WyioK t) = - <5^o(r, t), Yj/ q;0 (r, t)) Cff + 5 fc ( 5j , fe <z(r, t), Yj/ q;0 (f, *)> CH . 

Proof. Using (13.94)) and the approximation methods of Prop. 13. 2[ this follows im- 
mediately from the definition (|3.44|) of Lemma [3~T1 □ 

Observe that, unlike the Brillouin formula for the energy density (f3~89"l) . the 



approximation (|3.98H3T99| for the stress tensor does not involve the frequency dif- 
ferentiation of the susceptibility but simply its value at the given frequency u, a 
property discovered by L. Pitaevskii, [36], [221 §81], for the dielectric media. 
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4. Example: Maxwell's equations in an inhomogeneous TDD medium 

In this section, we apply the general construction developed above to the classical 
Maxwell equations in a material medium, 4, Section 1.1, Section 2.2], 



(4.1) 



I d t T> (r , t) = V x H (r, t) - 47rj cxt (f, t) 
\d t B (f,i) = -VxE(f,t), 

(4.2) V • B (r*, t) = 

in units with the speed of light c = 1. Here D, E, B, H are the electric induction, 
electric field, magnetic induction, and magnetic field respectively, which satisfy the 
following material relations 

(4.3) D (f, t) = E (f, t) + 4ttP (E; r, i) , B = H + 4ttM (H; r, t) , 

and j cxt is the external driving current. The one remaining Maxwell equation 

(4.4) V-D(r,t) =47rp cxt (r,t), 

with p cxt the external charge density, is automatically satisfied at all times provided 
it holds at a given time and that j cx t, Pcxt together satisfy the equation of continuity: 

(4.5) d tPext (r,t) + V-j e xt(r,t) = 0. 

We allow arbitrary external current j ex t, taking ()4.4j) as the definition of /9 ex t- 

We take the polarization P and magnetization M to be of the linear response 
form, [22, Chapter IX, Section 77], 

POO 

(4.6) 47rP(r,t) = (e(r) — 1) • E (r, t) + / xe (r, t) ■ E (r, t - r) dr, 

Jo 

poo 

(4.7) 47rM(r,t) = (/x(r) - 1) • H (f, t) + / X h (f, r) • H (f, t - r) dr, 

Jo 

with 

• e(f) and p(r) the static electric and magnetic permeability tensors, assumed 
to be real symmetric and uniformly bounded from above and below 

(4.8) e_l<e(f)<e+l, M-1<M^)<M+1, 

where 1 is the unit tensor, and e±, ji± > are constants. 

• Xe (r, t) and xh (r, r) the electric and magnetic susceptibility tensors, also 
real symmetric and satisfying a power dissipation condition, namely, 

(4.9) Im {Cxf (r, 0} > 0, C = w + if), r, > 0, 



Xf (r, C) 



/•OO 

/ e ict x F (r,t) di, F = E,H. 
Jo 



4.1. Hamiltonian structure of the field. We parameterize the field using the 
electric field D and the vector potential A as follows: 

(4.10) u = , n = -D, B = V A A 

in the phase space V = ^(R 3 , R 3 ) © ^(R 3 , M 3 ), where L| ff (R 3 ,R 3 ) denotes the 
space of vector fields X : R 3 — * R 3 , with inner product 

(4.11) (Xi,X 2 ) LL = / X 1 (r)-X 2 (r)d 3 r. 

4 " 4tt J r3 
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Thus for u as in (|4.10j) 

(4.12) || u f _ / jjn^i- -■ Ai n~\,.rF. 



^/ r >MI 2 + |a ( ^}. 



w 

We define the symplectic operator JonV with the matrix 

(4.13) j = (; - 1 

For the moment we do not impose a gauge condition on the vector potential A. 

We take II = — D to be the canonical momentum. The choice of sign puts the 
symplectic operator in the canonical form (|4.13|) . More important than the choice 
of sign is the choice of D as the momentum variable. This is essentially forced on us 
if we wish to use the formalism of i)3.1l and ^3.2[ since we should have a Lagrangian 
which does not depend on spatial derivatives of d t Q- This choice is also suggested 
by the coupling of an external current to the Maxwell equations and is in agreement 
with the standard Lagrangian density of relativistic field theory: L = F^F^" with 
Fy.v — d^Au — d u A^ where A v is the four vector potential [2~T] . It is a different 
convention, however, from that advocated by Sommerfeld [41] and adopted by us 
in our announcement of these results [5]. 

In a non-dispersive medium, with xe = Xh = 0, the material relations are 

(4.14) e(f) • E(r) = D(r), 

(4.15) \i (r) ■ H (r) = B (r) = V x A (r) . 

Identifying these equations with (jl.lObj) and recalling the classical expression for 
the electro- magnetic field energy in a static dielectric, i.e., 

(4.16) / {E(r) ■ e(r) ■ E(r) + H(r) • ■ H(r)} d 3 r, 
suggests parameterizing stress space with the vector 




(4.17) . - j 

Thus we take the stress space H = V 

The impedance K : V — > _ff is implicitly defined by (|4.15l BT4]) . since if u = f , 



(4.18) 
Thus 

(4.19) K = 
with 





K E (rO \ /-K E (r) 

K H (r)i v Vxj I K H (r)Vx 



(4.20) K E (r) = [vM^]" 1 and K H (r) = [vW]" 

well defined positive definite tensors by ((4 
The Hamiltonian is therefore 



(4.21) h = / h(r)dV, 

Jr 3 

with 

h(r) = i- {n(f) • [e(r)] ^ • n(f) + {V x A(r)} • [/i(r)] _1 • {V X A(r)}} 



3(5 
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The resulting equations of motion, expressed in the form (|1.10a[ [T.lObj) . are 
r IL(?,t)\ _ . (-V x {K H (r) •fk^i)}" 



(4 - 22) 9t \A(f,t)J 1 -KE(f)-f E (r-,i) 



(4.23) 



f E (f,*)\ / -K E (f)-U(f,t) 
f n (?,t)J K H (r>{VxA(f,()} 



(4.24) d t 



Eq. (|4.22p implies the two dynamical Maxwell equations once we take the curl of the 
second component and substitute II = — D, fE = Ke(^) _1 E and fn = Kn(r)~ 1 'H.. 
Similarly, the material relations (|4.141 I4.15[) follow from (j4.23p . The divergence 
condition (|4. 2[) is satisfied since V • {V x •} = 0, which also shows that dt V • D = 0, 
using (|4.22[) . Thus we may define the time independent external charge density 
Poxt(r) = • D(f, t), so that gU) and g3) hold (with j cxt = 0). 

When the system is driven by an external current j oxt , we replace (|4.22p by 

'll(r, t)\ /-V x {K H (r) • f H (f, t)}\ A^j oxt (r, i) N 
A(r,t)J { -K E (r)-t E (r,t) J + { 

The gauge freedom for A is related to the non-trivial kernel for the impedance: 

(4.25) kexK = {0} iJ grad 
with 

(4.26) ff grad = closure in L 2 (M 3 ; R 3 ) of {W (r) :ip(f)e C 1 (R 3 ) } . 

If (n(f,t), A(f,i)) is a solution to H4.22U4.23p or (|4TMll4T25j) . so is (U(f,t), A(r,t)+ 
X?ip(r)) for arbitrary (time independent) ip. This is essentially the same as the 
invariance under translation of the center of mass for mechanical systems (see £)5.2.2I 
below), with the significant difference that only the magnetic field B = V x A is 
directly observable, so we cannot detect the shift. Gauge fixing of A is implemented 
by the boundary condition A(r, t — —oo) = \7ijj G ker K att — —oo, with the choice 
of %l> having no effect on any quantity expressed in terms of B = V x A. Henceforth, 
we take if> = 0. 

4.2. Extended Hamiltonian for a TDD-Maxwell system. To relate the TDD 
dielectric medium to the general local TDD medium of we associate II = — D 
with the momentum p and A with the coordinate q. Respectively the electric field 
E is associated with / p and the magnetic field H is associated with / q . In a TDD 
medium, the material relations (|4.14l 14. 1 5[1 are replaced by 



poo 

e(r)-E(r) + / drx E (r, t) • E(r, t - r) = D(r) , 
Jo 

/■OO 

fi(r)-U(f,t) + / dTxn(r,T)-U(r,t-T) = VxA(f). 
Jo 



(4.27) 

(4.28) 
Defining 

(4-29) X it)m = ( ^ ' ' K !L> K ,-JfV 

puts the system exactly in the form (ll.llbp considered above 

n\ _ ft E \ , y 00 , , /f E 



(4.30) if I 2) {r, t) = I £ 1 (r, f) + jf X (r) I £1 (r, i - r)dr. 
Note that satisfies the power dissipation condition on H by (|4.9[) and (|4.8p . 
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The Hamiltonian for the resulting QHE is (after a permutation of coordinates): 
(4.31) H(U) = \(101,K\3) H = i(/C E U E ,/C E U E ) WE + i(/CHUH,/C H UH)H H , 



with 
(4.32) 



U = 



Vu H 

and extended impedance operator 



n(r) 

, U E = I 6> E (r,s) 



v<M?,s), 



A(r) 
U H = | H (r» 



'-K E (r) -T E N 
A V. = I 1 | , Kk = 
9, 



'Mr>Vx -T H N 



1 







o a 



The extended phase space and stress space are equal 

(4.34) H = V = V E ®V E , 
where 

(4.35) V E = Li r (E 8 ) R 3 )ei 2 (R,ii r (R 3 ,]R 3 ))ffii 2 (R,ii r (R 3 > K 3 )) 



The extended symplectic operator is 



p3 jm3\ 



Li 



p3 td)3 



(4.36) 



J = 



f° 








-1 













-1 














1 














1 
































-1 


V° 











1 


0/ 



The string coupling operators Tp are obtained from the susceptibilities as follows 

/oo 
^s)-<P F (f,s), F = E, H, 

with coupling functions 

(4.38) ? P (f,s) = i / cos(ws)? F (r»dw, F = E,H, 

where 



(4.39) ? f (t» = V2^K F (rO-Imx F (r,w)-K F (r), 



2w 



F = E,H 
F = E, 

F = H. 



The Hamiltonian (14.31[) is conveniently expressed as 



(4.40) 



H 



H(r)d 3 r, 
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with the density 



(4.41) H(r) = — <^ E(r) • e(r) ■ E(f) + H(r) • M (f) • H(r) 



|&<Mr, S )| 2 + \9v(r,s)\ 2 + \d s cf> H (r, s)\ 2 + |0 H (r, S )| 2 



ds 



Here E and H are related to the canonical variables II = — D, A, and ip E H by 

/oc 
te (r, s )-0 E (r, s )d s , 
-OO 



(4.43) H(r) = ^(r)- 1 • {V x A} (r) - [VM^]" 1 ■ / &(r, s) • H (r, s)d S 



As in (|3.1ip , the total energy is a sum of terms corresponding to the energy of the 
TDD system (the electromagnetic field) and the energy of the strings (the medium) , 
with no interaction term. This might be puzzling, however E and H, as defined in 
(|4.42[ 14~43|) . incorporate the interaction with the strings. 
The equations of motion for the extended system, 

(4.44) d t U = JX T /CU + R 

with R = col(— 47rj oxt , 0, 0, 0, 0, 0) the external current, may be written as follows: 
d t D(r, t) = V x H(r, t) - 47rj cxt (r, t) 
d t 9 E (r, s, t) = 9 2 E (r, s, t) + ?E (f, s) ■ ' E(r, *) 

d t (f> E (r,s,t) = 9 E (r,s,t) 
d t A(r,t) = -E{r,t) 



(4.45) 



G>i0 H (r, s, t) - 9 s V H (r, s, t) + ? H (r, s) • VKO • H(r, t) 
, dt<j>n{r,s,t) = 9n{r,s,t), 

with E and H given by [jO^Pg]) . 

4.3. Energy flux and stress tensor for the TDD Maxwell system. Study 
of the stress tensor in dispersive dielectric media has a rather long history, see 
[TT] . [38] . [E],El2] and references therein. In particular it is used to compute the 
ponderomotive and Abraham forces, [H], 0, EH §75, §81], P3 Section 2]. The first 
formula for the stress tensor was derived by L. Pitaevskii, [36] . [22l §81], [T5l Section 
3.2] for almost time harmonic fields in a transparent, i.e. lossless, medium. The 
formula was derived by applying thermodynamical methods and time averaging for 
a resonance circuit and a capacitor filled with the dielectric. Pitaevskii's formula 
is unexpectedly simple: one has to simply replace e and /i in the expression of 
the stress tensor for the case of non dispersive medium with respectively e (w) and 
/i (lu). This differs dramatically from the case of the energy density where one has 
to replace e and fi with nontrivial frequency derivatives [uje (oj)] and ^ [uifi (lu)]. 

In this section we treat the stress tensor for arbitrary fields — not necessarily 
almost-monochromatic — in TDD dielectric media, based on the formalism of 
We recover Pitaevskii's formula in the next section using Prop. 
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To make contact with the results of SJH we take coordinate variables q = A, 
4> = (0 E ,0 H ) and momentum variables p — II = — D, 9 = (0 e ,0h), with the 
spaces Vo , H® all equal to R 3 with the inner product 

(4.46) (v,w)v = (v,w) H o = (v,w)jjo 



47T 

The map K p is 

(4.47) K p (r)n = K E (r)-n, 

and K q = in (|3.2[) . asonly spatial derivatives of A appear in the Hamiltonian 
(|4~22]) . The maps Yi(r), i = 1, 2, 3, are as follows: 

(4.48) Yi(r)A = K H (r) • {e, x A} , 
with e, the unit vector in the i th coordinate direction, so that 

3 

(4.49) Y(r)-VA = ^Y^r^A = K H (r) • {V x A} . 

i=l 

The general representation (|3.9[) for the Lagrangian density of the extended 
system specializes in this case to 



(4.50) L(f,t) = i-||f E (r,t)| 2 + 2f E (r,t).|y" 



<; E (s,r)<l> E (s, r,t)ds 



/OO 
|5 t F ( S ,f,t)| 2 d S - \d s cf> F ( S ,r,t)\ 2 (Is 
...... . . -oo 



F— E.H 



with f p 


= f E and f q 


(4.51) 


f E (r,t) 


(4.52) 


f H (r,i) 



The vector potential A transforms as a vector under rotations, i.e., 
(4.53) g id A(r,t) = -e^^t) +e j A i (r,t), 

i.e., gi j = ei ® ej — ej (g) e^. The vectors f Ei n transform identically, g. 



p,q 

id 



Lemma 4.1. The system is homogeneous if and only if K E and F = E, H, are 

independent of f, and is isotropic if and only if they are scalars. 

Proof. This is obvious, except for the proof of (|3 . 32[) for Y given by (|4.48[) with 
scalar Kh, which is straightforward but tedious. □ 

Although we use the formalism of Sj3j we wish to express the resulting quantities 
using the usual electromagnetic field variables. We have already written the energy 
density (I4.41|l in this form. Using (|4.17|) . we identify the electric and magnetic 
fields and define, as for a non-dispersive medium, the magnetic induction B(r, t) — 
V x A(7% t) Thus, using the definition L = (dtQ,P) — H, we may express the 
Lagrangian density (|4.50|) as 

(4.54) L(r,i) = E(r,t)-r>(r,t) + — ]T / \0 F (s, r, t)\ 2 ds - H(r, t). 
n ''"f^e.h^- 00 
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The following theorem follows by elementary calculations 
Theorem 4.1. The energy flux vector for the extended Maxwell system (14. 50[) is 

(4.55) S(r,t) = i-E(f,t)xH(r,t). 

If the system is homogeneous and isotropic and perturbed by an external current 
j cxt (r, i) then the stress tensor corresponding to (|3.44[) is 

(4.56) T,j (r, t) = i- {H 4 B, + H^Bj + A. t [V x H] J + A,- [V x H] i } (r, t) 

+ ^hj | L (^ *) - ■ H(r, t) | + i {Aijexty - Aj ;j ext;i } , 

In view o/ (|4.54p and the equations of motion (|4.45p . T may 6e re-expressed as 

(4.57) T 4J (r, t) = -L {HjBj + H.B, - 2£ i;i H • B} (f, t) 

07T 

+ — {DjEj + DjEj - 2^-E ■ D} (r, i) 

87T 



| H (^*)-^ E l°° \0 F (s,f,t)\ 2 d s 
AijextyCr, i) + — ft {DjAj + DjA,} (r, i), 

07T 



where the Hamiltonian H zs given by (|4.41[) 

TTie Hamiltonian H and /_ |$f( s j t)| 2 ds, F = E,H, may be expressed as 
integrals over the history of the electro-magnetic field: 

(4.58) H(f;t) = - J {dA(r,t)+-E(r,t)-j cxt (r,t)}, 

J — OO 

and 



/OO 
\6 F { S ,r,t)\ 2 ds 
-OO 



= ap^l f I { aTXF ( fjil „i 2 ) +9TXF ( r - j 2t-i 1 -t 2 )}F(r,t 1 )-F(r,i 2 )diidi 2 , 

^ J — OO J — OO 

with a E (r) = s(r) =, a H (r) = n{r), and xf{t,t) = -xf^.-t), F = E,H. 

Remarks: i.) S is the familiar Poynting vector for the energy flux in a dielectric, 
ii.) The momentum density, by (|3.34| . is 

(4.60) Pi (r,t) = ~T> ■ 8iA(f,t) 



F=E,H 

with 

(4.61) * id (f,t) = ±-{D i A j -T> j A i }(f,t) 



1 f°° 

— 22 / e F -dMr,s,t)ds-d 3 <S> hJ (r,t), 

T7_T7 XJ OO 



1 f 00 

— 22 / { F;i0 F;J ■- F;j0F ;i } (^, S, i)ds. 

TP LT ** OO 



87F F=E,H 
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For a homogeneous system, the conservation law p.40p holds and we can express 
p in terms of the history of the electro-magnetic field by (j3.53|) , namely 

(4.62) Pi (f,t) = / {-d,T h3 (f,t') + [d t A} ■j oxt (r,t')}di'. 



The last term in ()4.57|) . dt {DjAj + D^A*} (r,t), is the time derivative of a 
symmetric tensor. We may drop it from the stress tensor provided we redefine the 
momentum density p i— > p + ^=dj{DiAj + DjAj}(r, t). Thus we may equally well 
take the following for the symmetric Maxwell stress tensor in a TDD dielectric, 



(4.63) T t ,j(r,t) = — {HjBj + H^-Bi + 6ij (H • piH. — 2H • B)} (r, t) 
+ — {D,Ej + DjE 4 + Sij (E • eE - 2E • D)} (f, t) 



1 f°° 

5i 'i-T- E / {l^Fl'-I^^Fl'K^^^ds + AUcxtu^t), 



F=E.H 



where we have started with (|4.57[) . dropped the last term on the r.h.s. and sub- 
stituted the expression (I4.41| for the energy density H (r , t) . The corresponding 
momentum density, i.e., r.h.s. of (|4.60[) + g^9j{DjAj + D_,Aj}(r, t), is 

(4.64) Pl (r,t) = -{DxB^f.tJ+ftrt^tjAjff.t) 

47T 

1 f°° 

F=E,H ,y_00 

where we have recalled that V • D = 47r / o cxt (by definition) and used the identity 

(4.65) D diA - D VA t = {D x B}< . 

In the above formulas for a TDD medium, D(r, t) ^ e(r,t)Ei(f,t) and B(r, t) 7^ 
fj,(r,t)H(r*,t). However, in the non-dissipative case when the susceptibilities in 
(|4.271l4~2"8|) vanish, xe = Xh = 0, the material relations reduce to D = eE and B = 
/iH and the last term in r.h.s. of (|4.63p involving F disappears. If, furthermore, 
there no are external charges or currents, p e xt = and j ex t;j = 0, then the formulas 
(14.631 |4TM|) turn into the familiar symmetric Maxwell stress tensor [5T} §33], 



1 

1 



(4.66) Tij(f,t) = — {EjDj + EjDj — SijE ■ D} (f, t) 



— {HjBj + B 3 H« - 8i.j (H • B)} (f, t), 

07T 



and momentum density 



(4.67) p(r,i) - i-DxB^t) . 

4.4. Brillouin formulas for the Maxwell energy density and stress tensor. 

As we have discussed in H3.31 one can derive rather simple formulas for the time av- 
eraged energy density and stress tensor produced by almost monochromatic waves. 
We refer to these formulas as Brillouin formulas, as it was Brillouin who introduced 
them for the TDD dielectrics, [221 §80]. In this section we present the specific form 
of these formulas for the electro-magnetic field in TDD dielectric media. 
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We remind the reader that the formulas are derived for almost harmonic waves as 
described in Q3.3\ We assume below, without comment, that we have a solution of 
Maxwell's equations with all fields in the form (|3.54[) describing a slowly modulated 
carrier wave of frequency uj; i.e., 

(4.68) E(r,i) = Re {e- iwt E n (r,t)} , H(f,i) = Re { C - ia "H (r, t)} , 

and similarly for D, B, where Eo, Ho, Do and Bo denote the slowly modulated 
amplitude of the wave. 

We start with the energy density H. Let us define the time averaged energy 
density with "no losses" (even if Im£p 7^ 0, F = E,H): 



(4.69) Hnl (r , t) 

= JL |E*(r,i) • ^[weftw)] ■ E (?, t) + K* (r, t) ■ ^N^)] ■ H ( 9 
where •* denotes complex conjugation and 

(4.70) e (f, w) = e(r) + RexE(r, uj), n (r, uj) = fi(r) + Rexn(r, uj). 



By p.89p . we see that Hnl is indeed the correct first order approximation to the 
time averaged energy density if the medium is lossless at uj: 



(4.71) H(r,t) = H NL (r,t) +0(6) (no losses), 

where <5 x i is the time scale over which the slowly varying amplitudes Eo, Ho, 
Do, Bo change noticeably, and it is assumed that 

(4.72) ImxE^w) = lmx}i(r,uj) = 0. 

In general the medium is absorbing at frequency u and (|4.72|) does not hold. As 
we have seen in i)3.3l there is in this case no simple expression for H. Instead, by 
(|3 . 72[) and (|3 . 88[) , we have an approximation for the time averaged power density, 

(4.73) 5jT(r,<) = -^-{E*(r,t)-ujlmxE(r,uj)--E (r,t) 



8tt 



Ho(r, t) ■ ujlrnxnir, uj) ■ Ho(f, t) 

Im [d t E* (f, t))} ■ jt-[colmx E (r, uj)} ■ E (f, t) 

lm[d t H*(r,t)} ■ ^-[ W Imx H (r»] • H (r,<) 



+d t H NL (r,i) + o{5). 

The last three terms on the r.h.s. involve time derivatives of the slowly varying 
amplitudes and are of order S. However, the first two terms, which are non-negative 
and describe steady dissipation to the medium, are of order 1 in general. Thus (see 

(4.74) m(r,t) = ^-{B* (r,t)-ujIm.x E (f,uj)--E (f,t) 

+ H*(r, t) ■ wlmxnCr, uj) ■ H (r, t)} + 0(6). 
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We now turn to the stress tensor T. The expression (|4.63|) was derived under 
the assumption of isotropy and homogeneity, so we suppose that 

(4.75) e(r)=£, (i(r) = fi, xv(r, w) = %e(w), xn{r, u) = xhM 

are position independent scalars. As in the general case treated in £!3.3l thc Brillouin 
formula for the time averaged Maxwell stress tensor is surprisingly simple. Using 
Prop. [3^41 to express the time average of the string Lagrangian, we see from (|4.63|) 
that the time averaged stress tensor (with no external current) is given by 

(4.76) J-(r, t) = ^jiic [E^D 0y + K, : D ::: , + ■ ■. 1 i. + H, :/ B, ; , (r, t) 



e(w) |E | 2 - 2Eq • D + |H | 2 - 2H^ • B [r,t) 



0(6), 



where e(u)) and h(uj) are the r independent versions of (|4.70|) . i.e., 

(4.77) e(w) = e + RexE(w), = /u + RexH^) • 

To simplify (|4.76|) even further, we use an approximation for the carrier wave am- 
plitudes Do and Bo, which is verified using the material relations (|4.27lfQ5]l . 

(4.78) D (f,t) - {e + XEH}E (f,t)+O(«5) 

(4.79) B (r,t) = { M + xh(«)} H (r,t) + 0{5). 
Thus, 



(4.80) TiAr,t) = 



16tt' 



e(w) 2RcE5. i E 0;j - <5 M - |E (r, t)| 2 (r,t) 



H(u)) 2ReHLH 0;j - 6 iti |H | 2 (r, t) + 0(6), 



Formula (|4T5D|) reproduces the Pitaevskii formula, [35], [H §81], [HI Section 
3.2] for the Maxwell stress tensor, derived in the references under the assumption 
of negligible losses at the carrier wave frequency u>. We note, however, that (|4.80|) 
is valid even if there are losses at ui! The main point of (|4.80[) is that in a TDD 
dielectric the Maxwell stress tensor has the same expression as in a lossless dielectric, 
with material constants incorporating the real part of the susceptibilities computed 
at the carrier wave frequency. This is in contrast to the energy (I4.7ip and power 
(|4.73jl densities, which involve frequency differentiation and, in the lossy case, the 
dissipative part of the susceptibilities. 

5. Precise formulation of the construction 

We now return to the general problem of constructing a QHE. The construction 
of Section [5] is correct, but is formal in two respects: 1) The Fourier transform of 
the susceptibility function may be an operator valued measure or distribution, in 
which case the point- wise limit (|2.1ip does not hold. 2) We have ignored domain 
questions. Specifically, we have not specified the domain of the extended impedance 
JC, nor have we shown that the dynamics of the extended system exists. 

Neither of these points poses a serious technical obstacle, and both are easily 
dealt with by established methods. We shall circumvent the first issue here by 
restricting ourselves to £( w ) defined point-wise almost everywhere. More general 
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susceptibilities could be handled by replacing the spaces L 2 (R, H ) with L 2 spaces 
with respect to an operator valued measure, via the Naimark construction [31] as 
in [7J. The reader familiar with the general theory can easily fill in the details. 

The second point is more essential, however, and will be dealt with carefully be- 
low. If the operator valued string coupling function c(s) is defined pointwise and is 
sufficiently integrable, then T is a bounded operator and this is relatively straight- 
forward. However, <;(s) may lack integrability or may be defined as a distribution, 
which might result in unbounded T. Thus we need to consider the definition of T 
and K, more carefully. 

We rely on some standard notions and results for operators on real Hilbert spaces, 
summarized in Appendix Kl We also use some notation defined there, in particular 

(5.1) H) — {closed densely defined operators from V — > H}, 

(5.2) B(V,H) = {bounded operators from V -► H}, 

with V and H real Hilbert spaces. We set C(V) = C(V, V) and B(V) = B(V, V). 

5.1. Hamiltonian evolution. The very first thing we require is that we can solve 
the evolution equations without dissipation. This is guaranteed by the following: 

Hamiltonian skew-adjoint condition (HSC): The symplectic 
operator J £ B(V) and impedance operator K £ C(V, H) are such 
that KJK T , defined on the domain 

(5.3) V(KJK T ) = {feV{K T ) : JK T feV(K)} , 
is skew-adjoint. 

Remark: Clearly KJK T is anti-symmetric on the domain (15. 3|) . To verify skew- 
adjointness we need to check that the domain is dense and the operator closed. 

The Hamiltonian skew-adjoint condition gives us a one-parameter group of or- 
thogonal transformations e tKJK on stress space H. We now show how to use this 
group to solve the non-dissipative initial value problem (see (|1.10a[ IT.10b|l l. 

(5.4) d t u(t) = JK T Ku(t), u(0)=u a . 

Because the generator JK T K may be unbounded, we do not try to solve (|5.4[) 
as such, but look for a finite energy weak solution u(t). That is, we seek a map 
1 1 ^ u(t) with u(0) = uq such that: i.) u(t) £ D(K) (finite energy), and ii.) 

(5.5) ^(u,Ju(t)) = -{Ku,Ku(t}), for any u G V{K) . 

In particular, we require the initial value to have finite energy, u G T>(K). 
To solve (|5.5I) . let f(t) — Ku(t) be the stress as in (|1 . 10h>[) and note that 

(5.6) d t (gj(t)) = -(KJK T gJ(t)), 
for any g £ V(KJK T ), so 

(5.7) f{t) = e tKJKT f(0) = e tKJKT Ku . 

That is, the stress is propagated by the orthogonal group e tKJK . The solution u 
may be obtained by integrating (JT7TJ) : 

(5.8) u(t) = u Q + JK T / /(t')dt', /(*) = e tKJKT Ku . 

Jo 
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The Hamiltonian skew-adjoint condition guarantees that J* f(t')dt' 6 V(K T ), since 



(5.9) KJK T f f(t')dt' = KJK T \ f e t,K JKT dt' 
is bounded by 2\\Ku \\ and V(KJK T ) C V(K T )E 

Theorem 5.1 (Constant energy evolution). Assume the Hamiltonian skew-adjoint 
condition holds for the pair J , K , and let u € D(K) be given. Then (|5.8|) is the 
unique finite energy weak solution u(t) € T^{K) to the initial value problem (|5 .4|) 
and the energy h(u(t)) = ^\\Ku(t)\\ 2 is a constant of the motion, h(u(t)) = h(uo) 
for all t. 

Proof. Taking u(t) as in (|5.8[) , clearly u(0) — uq and 

(5.10) J t ( v > Ju ^ = 
for v S V(K). Since 

(5.11) Ku{t) = Ku + [ d t ,f(t')dt' = fit), 



Ku = (e tKJK -l)Ku 







we see that uit) is a weak solution to 

Conservation of energy holds for any weak solution by (|5.11[) and the orthogonal- 
ity of e tKJK . To show uniqueness, it suffices to consider uq = 0, since the equations 
are linear. By conservation of energy, the solution u(t) € ker K at all times t, and 
hence by (|5.5p d t (v,u(t)} = for v in a dense set. Thus u{t) —u = {). □ 

The general solution to the driven Hamilton equation p,12[) is easily obtained 
by superposing solutions to the initial value problem (|5.4[) . noting that (|5.4[) is 
equivalent to (|1.12|) with p(t) = uo5(t) and u(t) = for t < 0. Thus, if we take a 
driving force p(t) € T)(K), the formal solution to (|1.12[) is given by 

(5.12) u(t-p) = f {p(t') + JK T f(t';p)}dt', 

J — OO 

(5.13) f(t-p) = f e^ KJKT Kp(t')dt', 



where lim f ^_ 00 w(t) = u_oo G kerif. Some assumption on p is necessary to guar- 
antee that (|5.121 I5.13P make sense. We shall require that the driving force was 
identically zero before some initial time, 

(5.14) p(t) = for t < t for some t e K, 
and that ||p(t)|| and are locally integrable 

(5.15) / {||/o(t)|| + ||ifp(t)||}dt < oo for any -oo < a < b < oo. 

J a 

In fact, (|5.14|) is overly strong as one only needs sufficient integrability at t — — oo. 
However, this assumption is convenient and not really restrictive from a physical 
standpoint. In any case, the r.h.s.'s of (|5.12i 15 . 13[) are well defined and u(t; p) is 
furthermore the unique weak solution to the driven Hamiltonian equations (|1.12p . 



It is key here that we have assumed that KJK T is closed on the domain specified in the HSC. 
If it were only closeable, we might not have T>(KJK T ) C T>(K T ). 
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Theorem 5.2 (Evolution under an external force). Assume the Hamiltonian skew- 
adjoint condition and suppose given an external force p{t) satisfying (|5.15l) and 
initial state at t = — oo. lt_oo G ker if. Then the unique weak solution to (|1.12|) 
with wk-linit^-oo u(t) = u_oo is given by (|5.12[ [5~13|) . 

Remark: A weak solution to ()1 . 12|1 is a function u(t) £ T)(K) satisfying 

(5.16) —{u,Ju{t)) = -{Ku,Ku{t)) + {u,Jp(t)), for any u € V{K) . 

When ker if is non-trivial, as for the electromagnetic field above, we have a gauge 
symmetry: if u(t) solves (|5.5[) or (|5.16[) then so does u[t) + vo with vo E ker if (for 
suitably modified initial condition in the case of (15. 5p ). With this symmetry comes 
a conserved quantity: the component of u in ker if J = J ker if. Indeed, by (|5.12D . 

(5.17) (u,u(t)) = (u,U-oo)+ [ (u, p{t'))dt' , ueJkerif. 



Thus, Pu(t), with P orthogonal projection onto J ker if, is constant, unless the 
driving force has a component in J ker if. As we have seen, in electromagnetism, 
translation by an element of ker if corresponds to a gauge transformation of the 
vector potential and projection onto J ker if singles out the electrostatic part of 
the electric field. 

5.2. Two examples. In the previous section we have presented a general abstract 
approach to quadratic Hamiltonian systems. Before turning to TDD systems and 
their extensions it may be useful to consider a couple of familiar examples viewed 
from the perspective of Thms. I5TT1 and I5T21 

5.2.1. A String. The vibrations of a Hilbert space valued string play a key role in 
the present paper, providing the dynamics of the auxiliary fields which give rise to 
the dispersion in the given TDD system. An H- valued string is also a good example 
of a system of the type analyzed in the previous section. The impedance operator 
and symplectic operator are 

(5.18) K = (I - J = (I "q 1 ) - onL 2 (R;H)®L 2 (R;H) , 

with V(K) = L 2 (R; H) V(d s ). This pair satisfies the HSC, with 

, T / 1\„ „ T ,/VT E rT\ 



(5.19) ifJif 1 = (J Q jds, on V(KJK l ) = V(d s )®V(d s ). 

The kernel of if is trivial. The orthogonal group e tKJK can be expressed in terms 
of the translation group e tds f(s) = f(s + t), i.e., 

tKJK^ 1 A lW , 1 ( 1 -l\„- fSs 



(5.20) e— = ye- + - ^ {)r 

Thus the solution to the initial value problem (|5.4[) is 

(5 - 2i) (Xm = (? ) (s) 

2 Jo I V °o + d s4>o ) y ' \ Oo- d s (j>o 1 
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Changing the s-derivatives into ^-derivatives and integrating gives 

(J)M-H( fc+ ^)" + " + jT*'U + o) 

with the two terms on the right hand side giving left and right traveling waves, 
respectively. 

5.2.2. Circular String. For the electromagnetic held, the zero modes corresponding 
to the kernel ker if of the impedance describe the gauge freedom of the magnetic 
potential and are not directly observable. In other systems these modes may be 
observable. An example of this type is provided by a circular string, with 

(5.23) K = (I °J , J = (J ~A , on L 2 (S 1 :H)®L 2 (S 1 ;H) , 
with V{K) = L 2 (S l ;H) ® V(d a ) and 

(5.24) KJK T = f? t)d a , onV{KJK T ) = V{d a )®V{d a ). 



1 / 

The solution (|5.8|l in this context is identical to (|5.22j) : 



(M,^{r + r)<-^r^u° + o)} 



(5.25) 



2 IV <fo / 2 _/„ \9 {a - f) 

The kernels of if and if J are the sets 

(5.26) ker if = i ( . J : 4>(a) = u= constant 



<jy{a) 



(5.27) ker if J = <j ( V Q ; ) : 0(a) = v = constant 



6(a) 


In this example, ker if _L ker if J and it is useful to write [ ,° ) as 

\<poJ 



(5.28) 



o/ V M 



\ / w \ , Wo 



wo , 



with f = J o (a)da, u — ^ J sl 0o(a)da, 6>o(a) = #o(a0 — u i 0o( a ) — u - In these 
coordinates, the solution (|5.25j) becomes the sum of clockwise and counterclockwise 
traveling waves superposed on uniform translation with velocity v. 



(5.29) \")(a,t) = 



v 

u + tv 
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where d^ 1 is a bounded left inverse to d a for functions / with J gl f(a)d9 — 0, e.g., 
d~ 1 f(a) = J " /(a')da'. The zero mode (0,u) T € kerif is the center of mass of 
the string and the conserved quantity (v, 0) T € Jker K is just it's velocity! 

5.3. Construction of a Hamiltonian extension. To construct the Hamiltonian 
extension we need two conditions on the susceptibility. The first is the 

Power dissipation condition (PDC): The susceptibility \{t) is 
a bounded symmetric operator on H for almost every t, 

(5.30) e~ et \\x(t)\\ S ^([O.oo)) , for every e > 0, 

and x(C) — a we M defined analytic function of C in the upper half 
plane by (|5.30p — satisfies the power dissipation condition ([1.14)1 . 

Eq. ([5.30P allows us to define the frequency domain susceptibility x(0> a ^ least 
for £ in the upper half plane. It is an analytic function taking values in the bounded 
operators on the complex Hubert space <CH. Since x{t) is rea l symmetric for a.e. t, 

(5-31) x(C) T = X(C), and x(0* = x(~C), 

where •* denotes complex conjugation. Thus x(0 — x(~C*) = X.(0* anc ^ 

(5.32) $(C) = ImCx(C) = | {CX(0 " C*X(0*} 

is a harmonic function of ^ in the upper half plane taking values in the bounded 
operators on the real Hilbert space H . By the PDC ([1.14[) , (/, $(£)/) is a positive 
harmonic function for any / e H . As such, a classic theorem of Herglotz states 
that there is a non- negative Borel measure fij onR such that 

(5-33) (fMOf) = - I |^^/(cM- 

The measure fi t in ([5.33P is the weak* limit as rj — * of the absolutely continuous 
measures (/, $(cj + vq)f)du). (In general there would be a linear term a /Im£ on the 
r.h.s. of (j5.33[) , however the PDC and dominated convergence imply that a/ = 
since ||$(C)|| = o(Im£) as Im£ — * oo.) 

The HSC and PDC together are sufficient for the existence of a quadratic Hamil- 
tonian extension. However, in general, this would require us to use the "L 2 space" 
with respect to the operator valued measure $(cj + i0)dw in place of L 2 (R, H) in 
the definition of the phase and stress spaces (jl.20ll.2ip . The construction of that 
space is not difficult (see [7] Appendix A] and references therein), but is abstract. 
The more concrete version of the extension, using the space L 2 (R, H), is available 
when the boundary measures are absolutely continuous. For this, we require the 

Susceptibility regularity condition (SRC): The measures fj,f 
appearing in ()5.33p are given by 

(5.34) p f (du) = (/,$( W + iO)/)d^, 
with $(cj + iO) the weak operator topology limit 

(5.35) $(w + i0) = WOT-lim^ *(w + i7?), 
assumed to exist and be bounded for almost every uel. 
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The power dissipation and susceptibility regularity conditions are general enough 
to cover many interesting examples. For example, if \ is of the form 

(5.36) X (t) = X°° + M*) 

with x°° a bounded strictly positive operator and ||/i(t)|| £ L 1 ([0,oo)), the condi- 
tions hold provided the continuous function 

(5.37) $(w + iO) = x°°+^ImMw) 

is everywhere positive semi-definite. The susceptibility could even grow as t — > oo: 
the function x(t) = t a for < a < 1/3 satisfies the PDC and SRC with 

(5-38) x(C) = T(a+1)^, 

and 

(5.39) $(w + iO) = T(a + 1) cos(mr/2)|w|- Q . 

The main application of the SRC is to the definition of the string coupling 
operator T : L 2 (R, H) — > H corresponding to 1T[) and (|2.2ip . A convenient way 
to organize this is as follows. By (|5.33[) and (|5.34[) . 

1 



(5.40) Sf(n) := — — + i0)/ , 

1 — IK 

defines an operator from H — > L 2 (R, CH), which is bounded since 

(5.41) \\Sf\\ 2 = 2 / ___</, $(« + i )/)d« = 27r(/,$(i)/). 

J-oo I 1 — lK \ 

Because 

(5.42) [£/(«)]* = Sf(-n) , 

the inverse Fourier transform of Sf is real (i.e., H- valued) almost everywhere. 
Thus, by the Plancherel theorem, 

l r R 

(5.43) Sf(s) := tf-limj^oo— / d Ke - 1KS S/( K ) 

defines a bounded map H — > i 2 (R, ) . We take the coupling operator T to be 

(5.44) T := S T (l-d s ) , 
on the domain T>{d s ) C L 2 (R, H), which is the set 

(5.45) V(ds) := {</> G L 2 (R, ff) : k$(k) £ L 2 (R, Cff)} . 

The coupling T may not be closed, nor even closeable, but it is densely defined with 
T(l — 9,,) _1 = S T everywhere defined and bounded. 

If the inverse Fourier transform of the coupling function 

(5.46) cr(s) = -!- / dKe- iK V 2$ (« + iO) 

exists in a suitable sense, for instance if || ^/2$(k + i0)|j £ L 2 (R) + L 1 (M), then 
(fl~T7| holds, i.e., 

/OO 
ds<,(s)<j)(s) . 
-OO 
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Even if (|5.47|) is purely formal, we may define the extended impedance 


1 



(5.48) 



K = 





K 





-S T 





f 











d s 



on the domain 
(5.49) 



V{JC) = V(K) 8 L 2 (R, H) © V(d a ) 



Theorem 5.3 (Skew-adjointness of the extended Hamiltonian). Assume the Hamil- 
tonian skew- adjoint, power dissipation, and susceptibility regularity conditions. Then 
the extended impedance operator IC, defined according to (|5.48|) on the domain 
(|5.49[) . and the extended symplectic operator J , defined by (|1.22|) . together satisfy 
the the Hamiltonian skew-adjoint condition. 



Proof. That K. is densely defined and closed is an easy consequence of (j5.48j) and 
the corresponding assertion for K and d s (on L 2 (M.,H)). So, the main point of 
the theorem is that /Cj77C T is skew-adjoint. To verify this, consider the domain of 
V{JC T ). Since 

(5.50) 



IC 1 











)• 


1 


1 










we have that 



(5.51) V{JC T ) = 





G V(K T )®L 2 (R, H)®V(d s 



In other words / G T>(K T ) and the linear combination Sf - 

<KJK T 
S 




G V(d s ). Thus, 












(5.52) 


K.JK. T = [0 


1 


•)( 




















)• 


1 












on the domain 



(5.53) V(K,JK? 




f G V(KJK T ), B, 



SfeV(d s ) 



The proof of skew-adjointness is now a simple exercise (see [3 Prop. 2.2]). 



□ 



5.4. Dissipative dynamics for the open subsystem. Combining Theorems 
15.11 15.21 and 15.31 we see that the Hamiltonian skew-adjoint, power dissipation, and 
susceptibility regularity conditions ensure that weak solutions to the extended sys- 
tem exist for suitable initial conditions or external forces. We now consider the 
relationship between these solutions and the initial dissipative system viewed as an 
open subsystem of the larger Hamiltonian system. 

We are interested in the relating the solution of and the solution of 

(5.54) d t U(t) = JIC T ICU(t)+I v p(t), 

with p : R — > V a driving force and ly the isometric injection of (|1.23|) . Recall that 

(5.55) JI V = WJ , and K.I V = I H K , 
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with Ih as in (11.23p . From the form of the impedance (15.48P it is clear that 

(5.56) ker/C = IykevK, 

which is to say that the equilibrium configurations of the extended system corre- 
spond to equilibrium configurations of the open subsystem. 

Suppose the extended system is driven out of an equilibrium state IyU-oo G ker JC 
at t = — oo, by an external force Iyp{t). The resulting trajectory is, by (|5.12|) . 

(5.57) U{t) = I v u_ 00 + j {l vP {t') + JIC T F(t')}dt', 

J — oo 

with the extended stress 

(5.58) F(t) = f e( t ~ t ')' c - 7 ' cT O y/9 (i')di' = f eW**? I H Kp(t')dt' . 
Thus the evolution of the open subsystem, described by u(t) — Iyll(t), is 

(5.59) u(t) = Uoo + f {p(t') + JK T f(t')}dt', 

J —OO 

where the kinematical stress f(t) is 

(5.60) /(*) = llF(t) = f lle^^I H Kp(t')dt' . 

The key point of the extension is that the pair (u(t), f(t)) given by (|5.59[[5T60)l is 
a solution to . More precisely, we interpret (| 1 . in the weak sense, namely 

(5.61a) d t (v, Ju(t)) = -(Kv, /(*)) + (v, Jp(t)), for any v G V(K), 

p OO 

(5.61b) (g,Ku(t)) = (g,f(t))+ (g, X (r)f(t - r))dr, for any g e H . 

Jo 

Theorem 5.4 (TDD dynamics for the open subsystem). Assume the Hamiltonian 
skew-adjoint, power dissipation, and susceptibility regularity conditions. Let S 
ker if and an external force 1 1— » p(t) S V satisfying (15. 14[) and (|5.15p be given. If 
u{t) = I$U(t) and f{t) = f&KU(t) are as in (15391) an d (OOP, then (u(t)J(t)) is 
the unique solution to (|1.11|) with wk-limt^_oo u{t) = Uoo- 

Proof. Let us first prove that (u, f) is a solution to (jl.lip . In fact (|5.61ap follows 
directly from (|5.59p . so it suffices to prove (|5.61b[) or, what is equivalent 

(5.62) (g, Ku(0) = (g, /(C)) + (3, 2(0/(0) , for any g G H and ImC > 0, 

where Ku and /, are the Fourier-Laplace transforms of Ku(t) and fit). Note that 
Ku and / are well defined since for t < to the external force p(t) = 0, by (I5.14p . 
and therefore Ku(t) = f(t) = 0. 

To simplify calculations, it is convenient to use the formalism of £12.11 in which 
the string coordinate is Fourier-transformed. This amounts to replacing K, by the 
operator appearing in (|2.28p . 



(5.63) 
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with k multiplication by the independent variable on L (R, CH) and T given by 



(5.64) T<t> = / v / 2<S>(n + iO)(j)(K)dK 



The domain of K, is 

(5.65) V(K) = V{K)®L 2 {K,H)®V{k), 
with 

(5.66) D(k) = : / K 2 ||0(K)|| 2 dK < 00 

I J —00 

Note that f is well defined for any cf S U(«) by (jQ3")l . 

Since the transformation /C — > /C results from the symplectic/orthogonal isomor- 
phism (|2.27|) , the relations (|5.59t [57^0]) hold with JC replaced by K. on the r.h.s., and 
a short calculation gives 

(5.67, /(C) - l{- 1 + lUj^ IJ ±—,„}K m , 

(5.68) Ku{() = 7- \-l + KJK T lJ I ^^- I h \kp(C), 

with p the Fourier-Laplace transform of the driving force p. Therefore 



1 



K 



(5.69) Ku(()-f(C) = 77 KJK'I^-I^JCJJC' ^-^——I H Kp{C) 



1 



KJK 1 + iC 



1 1 



Ti?(C) - / V2$(^+iO)i?(C;«;)dK, 



V^fiC %/2^iC 

where we have defined 

/ /(C) 

and noted that (see (|2.42[0 

(5.71) KJK T I%~ I^KJK T = (0 --^f 

Now, ([577D|l implies that 

K -K\(d(C-,K)\ _ 1 A/2$(« + iO)/(Q 



(5.72) , • ■ . . - , - — 



Thus 
(5.73) 



£(C;«A = J L_ fc V^ + iQ) /(C)\ 

£(C;«)/ v^^-C 2 ^V*(« + i0)/(Oy 
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The expressions for i9 and (p define square integrable functions of k, although 
v /2$(k + iO)/(C) may not be (see (|QD| ). Furthermore, by (|ol)9")) . 

f 

7T }-oo - 



(5.74) Ku(0-f(C) = -l ^-^$(« + iO)/(C)d K 



i r 00 r i f 



27r C J-oo l K ~ c -«- c 

i r°° i 



$(« + iO)/(C) d K 



-$(« + iO)/(C)d/s, 

""C ./-oo «~ C 

since 3>(k + iO) = <&(— k + iO). Finally, (]5.62[) follows because 

(5-75) -/ T<ff,*(K + iO)/)d« = C<0,x(O/>. 

for any /, g £ g 

To prove (|5.75p , it suffices, by polarization, to consider g = f. Then this identity 
is the well know representation of an analytic function with positive imaginary part 
and sufficient decay at infinity, see for instance [T] Section 59] or [Ml Section 32.3]. 
To see that it holds, note that the imaginary parts of the right and left hand sides 
agree by (|5.33p . Since the two sides are analytic functions of £, they can differ at 
most by a real constant, which must be zero as both sides vanish at £ = oo. 

It remains to prove uniqueness of the solution to (|I.11| . It suffices, by linearity, 
to consider Uoo — and p = 0. Then (|5.6Ia[) and (|5.62[) imply that 

(5.76) - (KJK T g, /(C)) = -i£(g, /(C)) - K(g, 2(0f) , 

for any g e V(KJK T ). Fix ( and let g = g x = A 2 |A + KJK T \- 2 f[Q. Then 
gx € T>(KJK T ) for A < oo and the l.h.s. of (|5.76[) vanishes since 

(5.77) (KJK T g x J(0) 

= x2 ( KJKT WJ^\ na WJ^\ Ti0 ) = °' 

since KJK T is skew-adjoint. Thus 

(5.78) = C(.9A,/(C)} + (.9A,CX(C)/(C)) 
which in the limit A — ► oo reduces to 

(5-79) CII/(C)II 2 - -(/(C),Cx(C)/(C)>. 

Taking imaginary parts, we see that this violates the power dissipation condition 
unless /(C) =0. As C was arbitrary / = and so / = and u = 0. □ 



Appendices 

Appendix A. Spectral theory of operators on a real Hilbert space 

In this section we review some known properties of linear operators on real 
Hilbert spaces. Our intention here is to set notation and to remind the reader of 
the differences and similarities between the real and complex cases by recalling the 
main results, without proofs. The material presented in this section is classical and 
well known to experts. Nonetheless, most of the standard textbooks focus on the 
complex case, making it difficult to point to a canonical source for the real case. 
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Recall that a real Hilbert space V is a vector space over the field of real numbers, 
complete in the norm topology induced by a symmetric scalar product (•, -) v . The 
norm of a vector u is ||u||y = y (u, u). We drop the subscript V from the norm 
and scalar product when it is clear from context which space we are discussing. We 
assume, without comment, that every Hilbert space considered is separable, i.e., 
has a countable dense subset. We denote the space of bounded linear maps from 
one real Hilbert space, V, to another, H, by B (V, H), and let B (V) = B (V, V). A 
linear operator K from V to a real Hilbert space H is a linear map K : T> (K) — ► H 
defined on a linear subspace T> (K) C V. The subspace T> (K) is the domain of K, 
and need not be closed. The range of K is the subspace, 

(A.l) K(K):={Ku:ueV(K)}CH, 

and also need not be closed. An operator K is densely defined if T> [K) is dense in 
V, and closed if T> (K) is a real Hilbert space, i.e., is complete, when endowed with 
the scalar product 

(A.2) (u, v) v{K) = (Ku, Kv) H + (u, v) v . 

An operator K is closeable if it has a closed extension, i.e., a closed linear operator 
K : V{K) -> H with V{K) D V (K) and Ku = Ku for u € V (K). The closure of 
a closeable operator fsT, denoted K, is the minimal closed extension, i.e., 

(A.3) V(K) = P| ^D(K) : K is a closed extension of k\ 

and Ku = Ku for any closed extension K. 

The transpose of a densely defined linear operator K from V to i? is a linear 
operator i£T T from H to V, defined as follows. Let V (if T ) be the set of vectors 
u' E H with the property that 

(A.4) \(u',Ku) H \<C u ,\\u\\ v , ueV(K). 

Since I? (K) is dense, the Riesz lemma — valid in a real Hilbert space, by the 
standard proof (e.g., see [39l Theorem II. 4]) — implies that for each v! £ T> (K T ) 
there is a unique vector K^u 1 such that 

(A.5) (u',Ku) H , = (K T u\u) , ueV(K). 

As it stands, V (K T ) might not be dense; indeed it might contain only the zero 
vector. However, as in the complex case, if K is closed, then V (K T ) is dense and 

Lemma A.l. A densely defined operator K from V to H , with V , H real Hilbert 
spaces, is closeable if and only if K T is densely defined, in which case K T = K 
and~K= (K T ) T . 

We denote the collection of closed densely defined operators from V to H by 
C {V, H). Thus - T :C(V,H)^C (H, V) and is an involution on C(V) = £ (V, V). 
Any bounded operator is closed and densely defined, i.e., B (V, H) C C (V, H). 

A complex Hilbert space V c is also a real Hilbert space, denoted here for emphasis 
by V^, under the inner product 

(A. 6) (u,v) v k = Re(u,v) Vi _ , 

where (•, -) v is the complex scalar product on V c - Note that || • \\ v k = || • \\v a , so 
these two spaces are identical as metric spaces. We denote the space of bounded, 
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respectively closed densely defined, complex linear operators from one complex 
Hilbert space, V c , to another, H c , by Be {V e , H c ), respectively Cc (V c , H c ). 

A complex linear operator K £ Cc(V c , H c ), considered as a map from Vf- to 
ii R , is obviously a real linear operator. However, not every real linear operator is 
complex linear. Indeed, a complex Hilbert space has a canonical symplectic operator 
Ju = iu (in the sense of (|1.5p ). such that a real linear operator K is complex linear 
if and only if it commutes with J. Conversely, given a real Hilbert space V and a 
symplectic operator J £ B (V), we may define a complex linear structure on V 

(A. 7) (a + ib) ■ u = au + bJu, 

and a complex inner product 

(A.8) (u,v) Vj :=(u,v) v -i(u,Jv) V) 

so that V becomes a complex Hilbert space, which we denote by Vj for emphasis. 

The complexification of a real Hilbert space V, denoted CV, is the complex 
Hilbert space, equal as a set to V © V, with multiplication by i given by 

(a.9) j = (; - 1 ) . 

An element u\ © 112, Uj € V, is denoted u\ + 1112, and we have 

(A. 10) (ui + iu 2 ,vi + iv 2 ) cv := (ui,«i) v - i {u 2 ,v\) v + i (u 1 ,v 2 ) v + {u 2 ,v 2 ) v 
(A. 11) (a + ib) ■ {u\ + 1U2) '■= au\ — bu 2 + i (bu± + au 2 ) . 

There is a natural operator of complex conjugation on the complexified Hilbert 
space <CV given by the block matrix 

(A.12) C = (J ° ± ) . 

The operator C is real linear, and, since CJ = —JC, is complex antilinear: Czu — 
z*Cu. We also use the notation Cu — u* , so 

(A. 13) (ui+iu 2 )* = U1—1U2, ui,u 2 eV. 

Any real linear operator A £ L (V, H) has a natural extension Ac G £; (CV, CH), 
called the complexification of A, namely 

(A.14) V(Ac) = {ui +iu 2 : Uj e V (A)} , Ac («i + m 2 ) = Am + iAu 2 - 

Not only is the extension Ac complex linear, it also commutes with conjugation 
AcC = CAc- Conversely any complex linear operator Ac £ Cc(CV,CH) which 
commutes with conjugation has the block matrix form 

(A.15) A C = 0) 

for a suitable real linear operator A £ £(V, H). Thus the category of real Hilbert 
spaces is equivalent to the category of complex Hilbert spaces furnished with dis- 
tinguished conjugation operators — self-adjoint complex anti-linear operators with 
C 2 = 1. From this viewpoint, real linear operators are complex linear operators 
which commute with conjugation. This leads to additional structure in the spectral 
theory, as we will see. 
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We show below, using the spectral theorem, that any symplectic operator J on 
V, has the canonical representation (|A.9jl by a block matrix for a suitable choice of 
real orthonormal basis on V. Associated to (|A.9|) is a real orthogonal decomposition 

(A. 16) V = Vi 8 Vi with dim Vi = dim V 2 , 

Given such a decomposition, we may interpret V\, Vi as "real" and "imaginary" 
subspaces of Vj respectively, with complex conjugation given by the block matrix 
(|A.12j) . (In the context of Hamiltonian systems V\ and Vi can be interpreted as 
the spaces of coordinate and momenta.) Note that the spaces in the decomposition 
(|A.16j) are not canonical: we cannot determine the real and imaginary subspaces 
uniquely from a complex structure. 

To develop a spectral theory for a real linear operator, it is generally necessary 
to work with its complexification. Indeed a real linear operator A G C (V) may 
have empty real spectrum. Thus we take the spectrum of A to be the spectrum 
of Ac, denoted a (A) = a (Ac), which is a non-empty closed subset of C for any 
A G C(V). (Recall that a (Ac) ={AeC | A c - A is not boundedly invertible }.) 

A curious phenomenon arises regarding the spectrum of a complex linear operator 
A 6 Cc (V c ). On the one hand we can define the spectrum as usual 

(A.17) a (A) = <r B(Vc) (A) = {A : A - A is not invertible in B (V c )} . 

On the other hand, we could forget that A is complex linear and consider V c as 
the real Hilbert space, Vj?, which we may then complexify to CV^, the set V c x V c 
with complex multiplication and inner product given by (|A.10| lA.llj) . Then the 
extension Ac of A has spectrum 

(A.18) a (A c ) = os(cv«) (A C ) = { A : A c - A is not invertible in B (CV C R ) } . 

The two sets (|A.17[ IA.18I) are not equal in general! Indeed, this is the case even 
for the map J = multiplication by i, since <Jb(v c ) (J) — {i} while a (Jc) = {ii}- 
However, this example already indicates the only difference that can occur: in 
general a (Ac) = a (A) U a (A)* . Indeed, 

Lemma A. 2 (Spectral Symmetry). Let A G C(V) with V a real Hilbert space. 
Then the spectrum of A, that is a (Ac), is a nonempty closed subset of the complex 
plane invariant under complex conjugation. 

Using the complexification, one can carry over to real linear operators various 
results from the spectral theory of complex linear operators. To state these results, 
we use the terminology that an isometry is a linear map T £ B (V, H) with ||Tu|| „ = 
\\u\\ v for all u G V, and an orthogonal map is an isometry T with 1Z(T) = H. 
(We reserve the term unitary for isomorphisms of complex Hilbert spaces.) An 
orthogonal map T is invertible with T _1 = T T . A partial isometry is a linear map 
T G B (V, H) which is an isometry when restricted to (kerT)" 1 , where 

(A.19) kerA = {u G V (A) : Au = 0} , 

for any linear operator A, and, for any S C V, 

(A.20) S ± = {u: (u, v) v =0 for all v G S} . 

In particular there is a spectral theorem for normal operators. Recall that an 
operator K G C(V) is called normal if K T K = KK T , i.e., D(K T K) = V(KK T ) 
and K T Ku = KK T u for any u in the common domain, and is called self-adjoint 
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if K T = K, i.e., V{K) = V (K T ) and Ku = K T u for all u G V(K). Clearly a 
self-adjoint operator is also normal. 

Theorem A.l (Spectral Theorem for Normal Operators). Let K G C(V) be a 
normal operator on a real Hilbert space V . Then there are i.) a a-finite measure 
space X,n, ii.) a fi-measure preserving involution $ : X — > X, i.e. $ o <&{x) = x, 
Hi.) a fi-measurable complex valued function k : X — > C, with 

(A. 21) k o $(a;) = k(x)* , fi a.e. x , 

and iv.) an isometry T : H — > L 2 (X,dfi;C), the complex Hilbert space of square 
integrable complex valued functions on X , such that 

(A.22) [TKu] (x) = k(x) [Tu] (x), \x a.e. x £ X, for any u G V (K) , 

and 

(A.23) ranT = {/ G L 2 (X, dfi; C) : / o = f(x)*, fi a.e. x) . 

Furthermore, {x : <&(x) — x} and {x : k(x) is real} differ only by a set of fj, 
measure zero. In particular, K is self-adjoint if and only if we may take $ to be 
the identity and k to be real valued. 

Thus the spectral theorem for self adjoint real operators is essentially the same 
as in the complex case. However, for a non-self-adjoint operator K, the involution 
<i> is non-trivial and represents additional structure compared to the complex case. 
This is related to the spectral symmetry described in Lemma IA.2| a strong version 
of which is eq. (|A.21|) . Associated to this symmetry is a natural partial symplectic 
operator J commuting with K, defined as follows. Consider the function a(x) = ±1 
if ±Imk(x) > and if lmk(x) = 0, and define a map J G B(V) by 

(A.24) [TJv]{x) := ia(x)[Tv]{x) . 

Because \a($(x)) = —ia(x), J is well defined. It is easy to see that Vo = ker J is 
the largest invariant subspace for K such that the restriction of K to this subspace 
is self adjoint. Thus = ker J 1 - is the largest invariant subspace on which K 
is "completely non self-adjoint." Furthermore, the restriction of J to Vq 1 is a 
symplectic operator. That is J is a partial symplectic operator, which is a bounded 
operator J such that J is a partial isometry and J T = — J . 

A special case of a normal operator is a symplectic operator J with J 2 = — 1 
and J T = — J. The spectral representation Theorem IA. II allows us to prove easily 
the canonical form (|A.9|) . Indeed, a symplectic operator being skew-adjoint has 
a spectral representation TJv(x) — j(x)Tv(x), with T : V — > L 2 (X, fi,C) an 
isometry. Since J 2 — —1, the function j takes values ±i. By (1A.23|) the following 
subspaces define an orthogonal decomposition of V: 

(A. 25) Vi = {v G V : Tv(x) is real for fi a.e. x} 

(A. 26) V% = {v G V : Tv(x) is imaginary for fi a.e. x} . 

Clearly J : Vi — > V% and J : V% — ► Vi. Furthermore, given a basis {ui,u 2 , ■ ■ .} for 
Vi, the sequence {Jui, Ju 2 , ■ ■ ■} is a basis for V2, and with respect to these bases 
the representation (|A.9|) holds. We have shown 

Lemma A. 3. Let J be a symplectic operator in a real Hilbert space V . Then there 
is a real orthogonal decomposition (|A.16j) in which the matrix of J has the canonical 
form (|A.9j) . 
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The functional calculus for a self adjoint operator associates to any real valued 
Borcl measurable function / : cr(K) -^la self adjoint operator / (K) with 

(A.27) \Tf{K)u]{x) = f(k(x))[Tu](x) , 

where the domain T>(f(K)) is the set of u such that the r.h.s. is square integrable. 
For a normal operator K, (IA.27P also defines a functional calculus, with f(K) a 
real operator if the map / : a(K) — > C satisfies f(z*) — f(z)*. 

The functional calculus allows us to define the square root and therefore the 
absolute value and polar decomposition, all of which work essentially as in the 
complex case. An operator K is called positive if K is self adjoint and 

(A.28) (u, Ku) v > 0, for all u £ V (K) . 

(For an operator in a real Hilbert space self-adjointness does not follow from (|A.28|) . 
so both conditions are necessary.) Given any K £ £(V,H) the operator K T K on 

V with domain V (K T K) = {u £ V (K) : Tu £ V (K T ) } is closed, densely defined 
and positive. (This is proved by the standard argument to obtain a self-adjoint 
operator from a semi-bounded quadratic form (see [391 Theorem VIII. 15]). A close 
reading of the proof shows that it works in a real Hilbert space.) Using the func- 
tional calculus we define 

(A.29) \K\ = VK T K, with \K\eC(V), V {\K\) = V (K ) . 

The polar decomposition of K is the factorization 
(A.30) K = T\K\, with T : V -> H a partial isometry, 

where T is uniquely defined if we require kerT = kerif, in which case 
(A.31) Tu = limK(\K\ + er 1 u. 

From the standpoint of the present work, the most important normal operators 
are the skew-adjoint operators. We call an operator K £ C(V) skew-adjoint if 
K = —K T . This implies that the function k given by the spectral theorem is purely 
imaginary, k{x) — ia{x)h{x), with o~{x) = ±1,0 and h(x) > 0. From the above 
discussion we see that the polar decomposition for a skew-adjoint operator is 

(A.32) K = J\K\ 

with J a partial symplectic operator commuting with K. Using J to put a complex 
structure on ran J, as in (|A.71 IA.8|) , we see that 

Theorem A. 2 (Canonical decomposition for a skew-adjoint operator). Let K £ 
C(V) be skew-adjoint. Then there is a direct sum decomposition V = Vq Vq with 
a complex structure on Vq (compatible with its real structure) such that K — 0©L4, 
with A a positive complex linear operator. 

The skew-adjoint operators are the generators of strongly continuous one pa- 
rameter groups of orthogonal operators, by the following analogue of the Stone- von 
Neumann theorem: 

Theorem A. 3 (Stone-von Neumann Theorem). A map t £ K — > Ot £ B(V), with 

V a real Hilbert space, is a strongly continuous, one parameter group of orthogonal 
operators if and only if 

(A.33) O t = e tA 

for a unique skew-adjoint operator A £ C(V). 
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Recall that Ot is strongly continuous if t i— ► Otv is a (norm) continuous map for 
any v £ V. 

Combining Theorems IA.2I and IA.3I we see that any strongly continuous one pa- 
rameter orthogonal group Ot decomposes as O t = 1 © XJ t with Ut a one parameter 
unitary group on a subspace. 

Appendix B. Sketch of the proof of Theorem 11.21 

In a remark following Theorem [L2] above, we briefly sketched a proof based on 
results from our previous work [Jj. To keep the present work self contained, in this 
section we give a more detailed sketch. 

To begin, note that it suffices to show 

Lemma B.l. The whole line extension (|1.37p of the friction function a e (t) can be 
written in the form 

(B.l) a e (t) = Te tLl T T , 

with L\ a skew adjoint operator on an auxilliary Hilbert space Hi and T : Hi — ► H 
a linear map. Furthermore, the minimal representation (jB.ljl is unique up to iso- 
morphism; i.e., if no proper subspace of Hi containing the range o/T T is invariant 
under L±, then L\ and T are unique up to transformation by an orthogonal map. 

Given the unique minimal representation (|B.1[) , we define C on H © Hi by 
(B.2) C = (_£ T Q . 

One may verify that this operator gives the desired extension, i.e., that (11. 4p holds. 
The uniqueness of the minimal extension follows from uniqueness in Lemma IB. II 

Strictly speaking the above argument works only when the 5 function contribu- 
tion 2a 00 S(t) to a e (t) vanishes, i.e., when x(0) = 0. If this term is non-zero, (|B.lj) 
holds with unbounded T, and some care must be taken in interpreting this relation 
(and also in the definition of C as in Thm. 15. 3p . It turns out that the required T is 
Li bounded, i.e., T : T>(Li) — ► H and 

(B.3) ||r/|| H < a\\L x f\\ Hx +b\\f\\ Hx 

for suitable finite a, b > 0. In place of (|B.lj) . we have a distributional limit: 

/oo 
(g(t),a c (t)f(t))dt 
-OO 

= ^fj9(thTe^(T*l) T f(t))dt, $i=0| + l) 1 

for all sufficiently smooth compactly supported H- valued maps / and g. The rep- 
resentation (|B.4j) and the corresponding uniqueness of the minimal representation 
are one direction of the following result related to an operator valued generalization 
of Bochner's Theorem [71 Theorem 7.1]: 

Theorem B.l. Let H be a real Hilbert space and let a e (t) — 2a oc 5(t) + a e {t), 
— oo < t < oo, with a e (t) a strongly continuous B(H)-valued function and a 
non-negative bounded operator on H . Then a e (t) is representable as (|B.4j) with L x 
a skew-adjoint operator on Hi andT an Li bounded map if an only ifa e (t) satisfies 
the power dissipation condition (|1 .36|) for every continuous H -valued function f(t) 
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with compact support. The operator T is bounded if and only if otoc = 0, in which 
case (|B.lj) holds. 

If the space H is minimal — in the sense that 

(B.5) {[Tf(Li)fv : / G C C (C) with f{z*) = /(*)* and v G h} 

is dense in Hi — then {Hi, L±,T} is uniquely determined up to orthogonal trans- 
formation. 

Sketch of the proof. The proof follows quite closely the proof of [71 Theorem 7.1], 
which deals with the complex case. We define the Hilbert space Hi to be a space 
of H value functions with the inner product given by 



(B.6) (<t>,i>) Hl = / / (<f>(t),a e (t-s)i)(s)) H dtds 



— oo •/ — oo 



The r.h.s. of (|B.6|) is non-negative by the power dissipation condition (|1.36| . but 
stritcly speaking only defines a semi- inner product. Thus we must mod out by the 
space of null vectors tp with = 0. 

The operator Li is then defined using the Stone- von Neumann Theorem IA.3I as 
the generator of the one parameter orthogonal group of translations 

(B.7) e sLl <t)(t) = <j>(t-s). 

Taking 

/>oo 

(B.8) Y<t> = / a e {~t)cf>(t)dt, 



it is not difficult to show that (|B.4[) holds. Note that, formally, 

(B.9) r T / = S(t)f. 

Uniqueness is shown as follows. Suppose we have a minimal representation such 
that (|B.5j) is dense in Hi, as can be obtained by restricting the constructed T and 
Lj to the closure of (|B.5j) . Given another representation 



(B.io) a e (t) = re tLi r T . 

with Li and Ti defined on Hi, one defines 

(B.ll) T[Tf(Li)ff = [f/(Z!)] T /, 

Then T extends to an isometric embedding T : Hi — > Hi, with 
(B.12) Tii = LT and TT = T. 

a 



Appendix C. Nonlinear systems with linear friction 

As mentioned in the introduction, the construction presented here extends to 
suitable non-linear systems. In this section we describe a few examples of this type. 

To keep the discussion simple, and to avoid the domain questions associated 
with unbounded functionals, we consider a finite dimensional mechanical system 
with V = R 2 " for some n, J a 2n x 2n matrix in the form (|1.5p . and Hamiltonian 

(C.l) h(it) = hi(u) + {Ku,Ku) H - 
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Here K is a linear map V — > H = R m , for some m, and hi is a C 1 function on V. 
Suppose the system evolves according to 

(C.2) d t u{t) = JK T f{t) + JVt,hi(u(t))+p(t) 

with dispersion in the definition of the kinematical stress /, 

poo 

(C.3) Ku(t) = f(t) + X (t- r)/(r)dr. 

Jo 

Here p(t) is a driving force and x( r ) is an m x m matrix valued susceptibility 
function which satisfies the power dissipation condition. 
Now consider the extended system with Hamiltonian 

(C.4) H(M,0) = -\\Ku-T<f>f H +hx(«) + 5 y {||*(*)|& + ||a s 0( S )|| 2 ff }d S , 

where = J_ <>(s)^>(s)ds with matrix valued string coupling function ? given in 
terms of x by (|2.22p . The resulting equations of motion, with driving force p, are 

(C.5) d t u(t) = JK T f{t) + JVhi(u) + p(t) 

(C.6) d t <f>(s,t) = 6{s,t) 

(C7) d t 6(s,t) = d 2 s cj ) ( S: t)+,( s ) T f(t) , 

where 

(C.8) /(<) = Ku{t) - Tcj>{t) . 

Note that given /, the equations for </> and 9 are identical to (|2.41 12. 3p . so the 
same driven wave equation (|2.6p for results. Thus the trajectory <f>, in terms of 
/, is given by (|2.7j) . resulting in the same equation (|2.8[) relating / and ifu. Since 
t is a solution to (|2.9p . we see that the 'w'-component of any solution to (jC.5IIC.8l) 
is a solution to (|C.2[ IC.3T) . The function hi played no roll in this argument. 

Of particular interest is a point particle subject to instantaneous friction. In this 
case u = (p,q) T , with p, q € M™ and J in the canonical form (|1.5p . and 

(C9) h(p,g) = V(g) + ^-(p,p) , 

with V(g) C 1 . To get instantaneous friction, we set the susceptibility 

(C.IO) X (r) = fa °q\ forallr>0, 

with 7 > the friction coefficient. Indeed, then 

iJ KJ V 1 J Wv(?(t))J 



(en) 0* (!)(*) 



with 

(C12) /(t) + 7 /°° /(* - r)dr = -Lp(i) . 

Jo V m 

Since /(t) = *Jmdtq(t) by (IC.lip . the "material relation" (|C 12[) implies 
(C13) d t q(t) +1 q{t) = -p(t). 

m 

Thus 

(C14) d 2 t q{t) = --VV(q(t))- 7 d t q(t), 

m 
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so q is subject to instantaneous linear damping. 

There is a Hamiltonian extension for this system with 



(C.15) H(p, 6,qA) = ^(p+V^mm) 2 +v(q) 

{(9( S )) 2 + (d s ^ S )) 2 }ds 



1 

+ 2 

and associated Lagrangian 

(C.16) L(q,d t q,<P,d t 4>) = -(d t q) 2 + v^^^O)) + - / (5 t 0( S )) 2 d S 

-V(«)-- J JdMfds 

Indeed, the Euler-Lagrange equations associated to (|C.16|) are 



(C.17) md 2 t q{t) = -W(g(t)) - y/2jTndt<l>(0,t) 

(C.18) 9?<KM) = d 5 2 0(M) + y/2^lS(s)d t q(t). 

The solution to (|C.18|) is easily seen to be <f>(s,t) = y/^^q(t — |s|), which when 
inserted into (|C.17|) gives (|C.14j) for q. 

The above ideas extend, with some care, to infinite dimensional systems. For 
instance, consider the non-linear wave equation 

(C.19) dfy(x, t) - d^(x, t) + V'Mx, t)) = >ydt%il>(x, t), 

with dissipation proportional to dtd 2 ip. By a formal extension of (|C.16j) this evo- 
lution could be seen as resulting from the reduction of the extended system 

(C.20) d t V(M) = d 2 V(M)- V'0/;(.t,<))- ^27 d t d x cf>(x,0,t) 

(C.21) d^(x,s,t) = d 2 <j)(x,s,t)- ^S(s)d t d x i>(x,t) 

with Lagrangian 

(C.22) Lfadt^hdtt) 

1 1 Z -00 1 

-|d t V(a:,*)| 2 + y/2^d t ip(x;t)d x <l>{x,0;t) + - J \d t <t>(x, s, t)| 2 ds| dx 

1 1 Z" 00 1 

-|d x #z,*)| 2 +V(#M)) + 2 J \d s <t>(x,s;t)\ 2 dsjdx. 

Of course, this begs the more difficult question of proving existence of solutions to 
the non- linear equations (|C 19[) and (|C.20p . 

Appendix D. Energy, momentum, and the stress tensor 

Consider a Lagrangian system in R d , described by a vector field Q(r) taking 
values in a Hilbert space Vq and with Lagrangian the integral of a density: 

(D.l) L(Q, d t Q) = [ L(Q{r),VQ{r),d t Q{r)-r)d d r. 

Given a field configuration Q : K d x [to, t{\ — ► Vq, the associated action is 

(D.2) A(\Q];ta,ti) = I ' / L(Q(r,t),VQ(r,t),d t Q(r,t);r)d d rdt, 
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The physical evolution Q(f, t) is a stationary point for the action A and thus satisfies 
the Euler-Lagrange equations 

(Recall the summation convention!) We use 5 to indicate partial differentiation of 
L to avoid confusion with dt and di, and write L(r, t), J^(r, i), ... as shorthand for 
L(Q(r, t), VQ(r, t), d t Q(r, t);r), ^(Q(f, t),VQ(f, t), d t Q(f, t); r), . . .. 

We have assumed that the Lagrangian density L does not depend explicitly 
on the time t. As a result, time translation is a symmetry of the system, and 
Noether's Theorem gives an expression for the energy density H(r, t), which is just 
the Hamiltonian, and the energy flux vector S(r) (see, for example, [2"Tl [5]): 

(D.4) H(r,<) = (d t Q(r,t),-^(r,t)^-L(r,t), 

and 



(D.:,i S,-(F./i = (d t Q(r,t),-^-(r,t) 



v„ 



When evaluated "on-shell," that is for Q satisfying (|D.3[) . these two quantities 
satisfy a local conservation law: 

(D.6) d t W{f,t) + diSi(r,t) = 0, 

implying that the integral of H 

(D.7) £ = I H(r,<)d d r 

is a conserved quantity, which we identify with the total energy of the system. 

Definition D.l. The system is homogeneous if L does not depend explicitly on r. 

A homogeneous system has an additional vector conserved quantitythe total 
wave momentum P. An initial expression for the wave momentum density is 

(D.8) pi(r,t) = (^iQ(r,t),^(r,t)J^ , 

with the associated canonical stress tensor 

(D.9) tij fat) - (w.iJ^K^ -SuLftt), 

see (2TJ El ED] ■ It is easy to show that 

(D.10) dtPftt) + djtijfrt) = 0, 

if the system is homogeneous, so 

(D.ll) P = / p(r,i)dV 

is a conserved quantity, which we identify with the total wave momentum. 
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As is well known [211 [3] , the conservation law and total momentum P are invari- 
ant under the following "gauge transformations:" 

(D.i2) Pi (f,t) = - OjQiAfy) 

(D.13) T itj (f,t) = T itj {f,t) + d t $ itj (f,t) + d k * itjtk {?,t), 

with 3>i j an arbitrary two tensor and ^i,j,k a 3-tensor antisymmetric in the later 
two indices: ^i,j, k = —^i,k,j- Indeed, for any such $ and W 

(D.14) / p(r,t)d d f = / p(r,t)d d r, 

and 

(D.15) ihp ■ ilj,., = dtp(f,t) + BjTij = 0. 

Thus the momentum density and stress tensor are not really uniquely defined. 
However, there is a well known way to fix these quantities, at least for an isotropic 
system. Here isotropy (see Defh. |P~2l below) expresses the invariance of the system 
under rotations and leads to another conserved current, the angular momentum 
density ^ , an anti-symmetric 2-tensor, which obeys the conservation law 

(D.16) d t m h3 (r) + d k F^ k (r) = 0, 

with F the angular momentum flux tensor. For isotropic systems we demand the 
following relations between the wave momentum and angular momentum densities 

(D.17) ™i,j{r,t) = r,pj(r,t) - r, Pi{n,t), 

and between the stress tensor T and the angular momentum flux tensor 

(D.18) F^r,*) = f;-T^(f,i)-^Ti )fc (r,t). 

For a homogeneous and isotropic system, the conservation laws (|D. 151 ID.16[) then 
imply that the stress tensor is symmetric, Tjj = Tjj, since 

(D.19) = dtmtj + d k F i>jtk 

= r-XdtPj + d k J hk ) - fj(dtPi + d k T hk ) +T jti -T id = T,,, T,. r 

The canonical stress tensor T is not symmetric in general, and thus is not the 
proper choice for a stress tensor related to the angular momentum flux tensor via 
(1D.18[) . For an isotropic system there are tensors $ and \& such that (|D.13|) is the 
symmetric stress tensor. 

To proceed we must define isotropy, and to do so must specify how the system 
transforms under rotations. Thus, we suppose given a representation T of the 
rotation group SO(cZ) by orthogonal operators on Vq. (See [5J 3D] for the basics 
of Lie groups and representation theory.) An arbitrary element of 50(d) can be 
expressed as e" with ui € so(d), the space of anti-symmetric d x d matrices, which 
is the Lie algebra of SO(d). Thus, the representation T : SO(<i) — > B(Vo) can be 
written in terms of a corresponding representation v of so(d) 

(D.20) T(e w ) = c^H 

The matrices ey € so(g0, 

(D.21) e i,j;a,f3 = 0~i,aO~j,f3 ~ 8j, a 8i,p, 
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satisfy eij — and the collection {e;.j : 1 < i < j < d} is a basis for so(d). An 
arbitrary element uj € so(d) can thus be written 

(D.22) ui = —uji^eij (summation convention). 

Thus the representation v is determined by the skew-adjoint operators 
(D.23) G id = v(e itj ), 

which evidently satisfy Gi.j = —Gj^ and 

(D.24) [Gi,j,G kil ] = v([e lJ: e k j}) = -S itk G jtl + 8 jt k G it t + S it iG jtk - SjjG iik . 

We assume that v is a representation by bounded operators, Gjj € B(Vb) for every 
The representative of a generic element ui € so(d) is therefore 

(D.25) v(uj) = -ui^Gij, 

and of a rotation e w 6 SO(<i) 

(D.26) T(e") = e 5 Wi ^ G M. 

A global rotation of the coordinate system about a given point r involves a 
transformation of r 

(D.27) r k f" = r° + e" • (f-r°), 

and of the field Q: 

(D.28) Q(r) ^ Q"{r) = e^- G -Q(r-"). 

Such rotations form a representation of SO(d) in L 2 (R d , Vq) with generators 
(D.29) L r ;°Q(f) - G hJ Q(r) - (n - f?) d,Q(r) + (r, - r?)diQ(r), 

so may be written 

(D.30) Q w (r) = [e^-<°Q](r). 

Definition D.2. The system is isotropic at a point r if the Lagrangian density 
at the given point r is invariant under the transformations (|D.30P : 

(D.31) L(e^'^^Q{r,t),d t e^'^ L ^Q(r,t),ye^'^ L ^Q(r,t)) 

= L(Q(f,t),d t Q(r,t).yQ{f,t)), for every uj eso(d). 

We say the system is isotropic if it is isotropic at every point. 

It is convenient to express (ID.31|) in infinitesimal form by differentiating the l.h.s. 
at ui = 0. To this end, note that 

(D.32) e^ L lsQ{r,t) = e^ G *->Q(r,t), 

since the remaining terms vanish at the origin r of the rotation. Similarly, 

(D.33) d t ei"^ L liQ{r,t) = e^ G ^W t Q{r,t). 

However, 

(D.34) d k Ll d Q(r, t) = G^d k Q(f, t) + 5 jtk diQ{r, t) - 8 hk d Q{r, i), 

because 

(D.35) [dk,^] = 8j, k di - Si^dj ^ 0. 
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Lemma D.l. The system is isotropic at r if and only if for any Q 

I S\- 
(D.36) = TiA^-T^t) + (G tJ Q(r,t),—(f,t) 



SQ 



with T given by (|D.9j) . If Q satisfies the Euler- Lagrange equations (|D.3|) then 
(D.37) = f itj (f,t)-fj,i(r,t) 



+ ft (GtjQVt), ^(r,t)) + dk (GuQftt), ^(r,t) 

Proof. Eq. (|D~36f is a consequence of the definition (|LX3T|) and (|1T32hD\341) . Eq. 
(|D.37[) follows from rewriting the third term on the r.h.s. of (|D.36[) using the Euler- 
Lagrange equations (|D.3[) and combining terms with the Leibniz rule. □ 

Following [5J Section III. 4] we define the wave momentum density and stress 
tensor by gauge transformations (|D.12I [D.13[) of p and T, with 

(D.38) <Mr,i) = l(GuQ(W>^ftt)) v > 

and 



(D.39) * idtk (?,t) = \\( G idQ(r,t),-^(m 



V 



Vq ) 

Note that ^i,j,k is anti-symmetric under interchange of j and fc. 

Theorem D.l. Let the stress tensor T be defined 

(D.40) Tij(f,t) := f id (r,t)+d t $ij(r,t)+d k V idtk (r,t), 

with <!>, ^ given by (|D.381 [D.39[) . If the system is isotropic, then T is symmetric. 
If Q satisfies the Euler- Lagrange equations and the system is homogeneous, then 

(D.41) d tPi (r,t) + (>,T,., I r. I) = 0, 

with the wave momentum density 

(D.42) Pi (r,t) = Pi (r,t) - dj^ji^t). 

If the system is homogeneous and isotropic, then the local angular momentum con- 
servation law (|D.16|) holds with m and F defined by (|D.17[) and (|D.18[) respectively. 

Proof. The only point not established in the above discussion is the symmetry of 
T. This however follows from Lemma ID. II since 

(D.43) Tij(r,t)-T Jyi (r,t) = r.h.s. of flE37|). □ 

To close, we consider how the continuity equations (|D.6UD.4"T]) are modified by a 
driving force R(r, t) € Vq such that Q satisfies the driven Euler-Lagrange equation 

< D ' 44 > *^«0 + *^('.0-|g(»W-««'>. 
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which is simply the Euler-Lagrange equation for the time dependent Lagrangian 
density L + (Q(r, t), R(r, t)}. Since R breaks time translation invariance and homo- 
geneity, energy and momentum are no longer conserved. However, we have 

Theorem D.2. If Q satisfies the driven Euler-Lagrange equations (|D.44I) then 

(D.45) dtHftQ + djSjftt) = {Q{r,t),R(r,t)) Va . 

If, furthermore, the system is homogeneous then 

(D.46) d tPi (r,t) + djTiji^t) = (diQ(r,t),R(r,t)) Vo . 

Proof. This is a straightforward calculation. □ 
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